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ABSTRACT 

We extend our investigation of magnetic field evolution in three-dimensional flows driven by the stationary 
accretion shock instability (SASI) with a suite of higher-resolution idealized models of the post-bounce core- 
collapse supernova environment. Our magnetohydrodynamic simulations vary in initial magnetic field strength, 
rotation rate, and grid resolution. Vigorous SASI-driven turbulence inside the shock amplifies magnetic fields 
exponentially; but while the amplified fields reduce the kinetic energy of small-scale flows, they do not seem 
to affect the global shock dynamics. The growth rate and final magnitude of the magnetic energy are very 
sensitive to grid resolution, and both are underestimated by the simulations. Nevertheless our simulations 
suggest that neutron star magnetic fields exceeding 10 14 G can result from dynamics driven by the SASI, even 
for non-rotating progenitors. 

Subject headings: magnetohydrodynamics (MHD) — methods: numerical — physical processes: turbulence 
— stars: magnetic field — supernovae: general 
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1. INTRODUCTION 

Not long after the discovery of pulsars — who se charac- 
teristic signal was linked to magnetic fields ( |Hewish et al.| 
1968) — the potential role of magnetic fields in the core- 
collapse supernova (CCSN) explosion mechanism began to 



-.41 



be investigated (e.g., |LeBlanc & Wilson 1 1970} [B isnovatyi- 
|Kogan et al. "[ [T976; Mei er et al. ||l<)76[|Sym balisty 1984). 
In principle, a differentially rotating proto-neutron star (PNS) 
could both amplify magnetic fields and serve as an energy 
reservoir available to be tapped by those fields, giving rise to 
magnetically powered explosions. An early conclusion, how- 
ever, was that both unrealistically rapid rotation and unreal- 
istically strong magnetic fields would be needed at the pre- 
collapse stage for magnetic fields to play a principal role in 



the explos ion dynamics (LeBl anc & Wilson||1970 Symbal- 
|isty |1984| >. 

In more recent years interest in strong magnetic fields has 
returned in connection with a number of observables re- 
lated to core-collapse supernovae, including asymmetries in 
the explosion ejecta (Wheeler et al. 2002), natal neutron 
star kick velocities (Lai & Qian 1998), and especially the 
high-energy electromagnetic activity connected to some neu- 
tron stars known as magnetars, or Anomalous X-ray Pulsars 
(AXPs) and Soft Gamma Repeaters (SGRs) (e.g., IDuncan 



& Tho mpson 

2005 
h 



Hu 



al. 



1992||Thompson & Duncan | 2001~] [Hurley et 
Woods & Thompson 2006). AXPs and SGRs 



10 35 erg s _1 , with sporadic outbursts releasing up to 10 4i erg 
per event. Gamma-ray outbursts from SGRs are even more 
energetic, an extreme example being the giant flare from SGR 
1806-20, which releas ed an estimated 10 46 erg over 380 s 
([Hurley et al. ||2005|l. Furthermore, AXPs and SGRs are 



neutron stars characterized by relatively long rotation periods 
(P > 1 s) and high spin-down rates (P > 10~ 12 ss _1 ) (e.g. 



Lorimer & Kramer 2005 ). As their rotational energy cannot 
account for the electromagnetic emission, and because of the 
strong magnetic torques implied by high spin-down rates, they 
are believed to be young neutron stars powered by dissipa- 
tion of extremely strong surface magnetic fields (10 14 -10 15 G, 



Duncan & Th ompson |1996 Thompson & Duncan 2001 ). 

On the theoretical side, trie discovery of the magneto- 
rotational instability (MRI) by |Balbus & Hawley |(|1991|l and 
its application to CCSNe (initiated by Akiyama et al. 2003) 
relaxed the requirement of strong pre-collapse -/^-fields, re- 
newing interest in magnetic fields as a possible key ingredient 
in the explosion mechanism of some supernovae (i.e., those 



from rapidly rotating progenito r cores; e.g., [Wheeler et al. 
2002 ; Obergaulinger et al. 200 5)|Moiseenkoetal. 2006 ; Bur- 
Irowsetal. 2007 ; Takiwaki et al. 2009J. (The MRI results 
in exponential growth of the magnetic energy on the rotation 
timescale.) However, the rotational energy falls off quadrati- 
cally with increasing rotation period, and is about 5 x 10 49 erg 
for a 20 ms period PNS — much less than the characteristic 

51 



are characterized by quiescent X-ray luminosities as high as CCSN explosion energy of - 10 5i erg = 1 Bethe (B). Thus 

any magneto-rotationally driven supernovae likely would be 
peculiar events, since magnetic progenitor cores tend to rotate 
slowly at the pre-collapse stage ( Heger et al. |2005} . 

Leaving aside the explosion mechanism, the relationship 
between the formation of neutron star magnetic fields and CC- 
SNe is still an open and interesting question, particul arly in 
the ca s e of magnetars (AXPs and SG Rs) (Lorime r & Kramer| 
2005 ). Thompson & Duncan ( 1993 i argued that such strong 
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fields must be generated during the neutrino cooling epoch 
after the collapse of the progenitor's iron core, and possibly 
before the explosion is initiated (~ 1 s after core collapse). 
Their model remains one of the prevailing theories for mag- 
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netar formation, and includes a convective a — Q, dynamo, 
which operates when the rotation period is comparable to the 
turnover time of entropy-driven convection (< 3 ms) near the 
surface of the PNS. The rapid turnover time may suggest that 
magnetars are formed in the magnetically-driven explosion of 
collapsed, rapidly rotating progenitors, whose remnant is spun 



down by MHD processes at later times. Bonannoetal. (2003 
2005) found that neutron finger instabilities (e.g., IBruenn & 



Dineva 1996 ) may also result in dynamo action in PNSs with 
rotation periods as long as 1 s. In this scenario, the formation 
of neutron star magnetic fields may be slow (compared to the 
explosion time scale), and their creation is not necessarily tied 
to dynamics in the supernova explosion. The MRI may also 
operate near the surface of the PNS, and contribute to neutron 
star magnetization. 

The lack of sufficient rotational energy in magnetized pre- 
collapse progenitor cores, as predicted by stellar evolution 
models (Heger et al. ||2005| ), has sparked some recent inter- 
est in MHD processes in non-rotating CCSN environments 



(Endeve et al. 2010 Guilet et al. 201 1 Obergaulinger & 



anka 201 l| i. These studies investigate field amplification 
mechanisms and the possible role of amplified 5-fields on the 
dynamics of slowly or non-rotating collapsed progenitors, in 
wh ich rotational MHD pr ocesses are insignificant. In particu- 
lar, Endeve et al. (2010 hereafter Paper 1 1 studied magnetic 
field amplification by the stationary accretion shock instabil- 
ity (SASI, Blo ndin et al.||20"03"l >. The SASI is central to the 
theory of CCSNe: recent simulations lead to the conclusion 
that it lik ely plays an imp ortant rol e in neutri no-powered ex- 
plosions (Bruennet al. 2006 ; Buras et al. 2006 ; Mezzacappa 



itlT1|2b07| IScheck et al. ||200tj| |Marek & Janka ||2009l 



Suwa et al. |2010 Miiller et al. |2012 1, and may also explain 
certain observables of pulsars , including their proper motion 
( IScheck et al. |2004[ > and spin ( |Blondin & Mezzacappa |2007| l. 
Thus, magnetic fields may be an important part of a supernova 
model if the S ASI is found to be sensitive to their presence. 
In Paper I we adopted the idealized model of |Blondin et 



al. (2003) and Blondin & Mezzacappa (2007 1, and added a 
weak radial (split monopole) magnetic field. We presented 
results from 2D (axisymmetric) and 3D MHD simulations of 
the SASI, and found that SASTdriven flows may result in sig- 
nificant magnetic field amplification. Magnetic field evolu- 
tion in axisymmetric simulations was found to be geomet- 
rically constrai ned. Moreover, the non -axisymmetric spiral 
SASI mode (Blondin & Mezzacappa 2007) dominates the 
post-shock flows in 3D simulations at late times. The non- 
linear evolution of the spiral mode drives vigorous turbulence 
below the shock, which results in exp onential am plification 
of B-fields due to "stretching" (e.g., IQtt ||1998[ ), and the 
magnetic energy becomes concentrated in intense, intermit- 
tent magnetic flux ropes. We presented results from models 
with non-rotating and weakly rotating initial conditions, and 
weak (10 10 G) and stronger (10 12 G) initial magnetic fields. 
The magnetic fields were not found to reach dynamically sig- 
nificant levels (i.e., components of the Maxwell stress ten- 
sor did not contribute significantly to the total stress), and 
hence no impact of magnetic fields on local or global dy- 
namics was demonstrated. However, we found that SASI- 
induced turbulent magnetic field amplification is very sensi- 
tive to the spatial resolution adopted in the numerical simula- 
tions. Most of the 3D models presented in Paper I were per- 
formed at "medium" spatial resolution (grid cells with sides 
Al = 1.56 km), while one model was performed with "high" 
spatial resolution (Al = 1.17 km). The thickness of magnetic 



flux ropes was found to decrease in proportion to Al. We 
did not observe convergence of i?-field amplification with in- 
creasing spatial resolution. Nevertheless, the simulations im- 
plied neutron star magnetization as a result of SASTinduced 
magnetic field amplification. 

This paper continues and extends the investigations initi- 
ated in Paper I It improves on our previous study in several 
important ways, including (1) coverage of a larger parame- 
ter space, (2) higher spatial resolution (up to 1280 3 zones), 
and (3) computation of kinetic and magnetic energy spectra. 
With the new set of simulations we investigate in some de- 
tail the nature of S ASI-driven turbulence, and the growth and 
impact of magnetic fields during operation of the SASI. We 
investigate the saturation level of magnetic energy in our sim- 
ulations, the (kinetic) energy reservoir available for magnetic 
field amplification, and the factors determining the magnetic 
energy growth rate. We also consider as in Paper I the im- 
pact of amplified magnetic fields on global shock dynamics, 
in particular any impact they may have on the SASI. Finally, 
we attempt to quantify the levels of neutron star magnetiza- 
tion that may be expected from SASI dynamics. 

We find that the SASI-driven turbulence shares several sim- 
ilarities with non-helical turbulence (e.g., Brandenburg et al.| 



1996; Haugen et al. 2004 1, and results in an efficient 



small-scale dynamo. Magnetic fields grow exponentially in 
the turbulent flows driven by the SASI as long as the "kine- 
matic regime" remains valid. The kinematic regime ends 
when the magnetic energy becomes comparable (locally) to 
the kinetic energy of the turbulent flows — the magnetic en- 
ergy source. From the computed energy spectra we estimate 
the "turbulent" kinetic energy, E 1 ™, available for magnetic 
field amplification, and, in our idealized model, E£ consti- 
tutes about 10% of the total kinetic energy below the shock 
(E ki „ ~ 5x 10 49 erg). The total magnetic energy saturates 
at about E mi , e ~ 5x 10 47 erg. The presence of amplified 
magnetic fields results in less kinetic energy on small spatial 
scales, but we find no impact of magnetic fields on global 
shock dynamics, which is consistent with considerations of 
the energetics. However, magnetic field evolution remains 
sensitive to numerical resolution, and magnetic fields are sub- 
ject to significant numerical dissipation during the saturated 
state, and our ability to quantify fully the impact of magnetic 
fields in a more realistic situation is therefore limited. The 
magnetic energy growth time decreases with increasing reso- 
lution, and, based on the turnover time of the SASI-driven tur- 
bulence, is estimated to be a few milliseconds. We argue that 
the MHD processes studied in this paper may contribute sig- 
nificantly to strong, small-scale neutron star magnetic fields, 
and provide a connection between the magnetic fields of neu- 
tron stars at birth and supernova dynamics. The saturation 
energies may be sufficient to power flaring activity of AXPs, 
and possibly SGRs. Moreover, their formation does not re- 
quire progenitor rotation. 

2. SETUP OF NUMERICAL SIMULATIONS 

We employ the same numerical code and three-dimensional 
initial conditions we used in |Paper I| which fol l ow clo sely 
the adiabatic se tup desc ribed in |Blondin et aL| ( |2003j l and 
|Blondi n & Mezzacappa (2007J: a stationary accretion shock 
is placed at a radius r = R Sh = 200 km, and a highly 
supersonic flow is nearly free-falling towards the shock for 
r > R sh with e ti „ + e mv w 0. Between the shock and the PNS 
the flow settles subsonically — obeying the Bernoulli equation 



P 



— and is nearly in hydrostatic equi- 
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Figure 1. Plot of the initial condition for the non-rotating weak-field model 
with So = 10 10 G (BIOLO.OEOO): internal energy density (e in „ solid line), 
magnitude of gravitational potential energy density (|e grav |, dash-dot line), 
kinetic energy density (e tirl , dotted line), and magnetic energy density (e mag , 
dashed line) versus radial distance from the center of the PNS. The surface of 
the PNS (our inner boundary) is fixed at r = R ¥NS = 40 km and the shock 
is initially located at r = R S h = 200 km. Inside the shock |eg, av |, eim and 
e mi , g follow roughly the same power-law (oc r~ 4 ), while e t \ n <x r . The 
flow is in steady state free-fall outside /? sh , with e kin and |e gr „| proportional 
to r~ 2 ' 5 , and e in , <x r -2 . e milg has been multiplied by 10 6 to become visible 
on the plot. (The dashed line is also identical to e mag in the strong-field model 
(B13L0.0E00), cf. Table[T]) 



librium. Matter is allowed to flow through an inner boundary 
placed at r = R PNS — 40 km. The mass density and pres- 
sure just inside R ms are determined from values just outside 
i? PNS using power-law extrapolations: p ex r~ 3 and P ex r~ 4 , 
respectively (a procedure that proved necessary in order to 
maintain the steady state of the unperturbed initial condition). 

Figure [T] displays the initial configuration of a spherically 
symmetric, non-rotating stationary accretion shock with a 
weak radial magnetic field (B = 1 x 10 10 G; the initial mag- 
netic fields in our simulations are discussed in further detail 
below). We plot internal energy density e in , = P/(-y — 1), 
kinetic energy density e ki „ = pu • u/2, magnetic energy den- 
sity e mag = B • B/(2/io), and the magnitude of the gravita- 
tional potential energy density e grav = p<I> versus radial dis- 
tance from the center of the star. Here p, u, P, B, and $ are 
the mass density, fluid velocity, fluid pressure, magnetic flux 
density (magnetic field), and gravitational potential, respec- 
tively. The vacuum permeability is denoted p,Q. We adopt the 
ideal gas equation of state, with the ratio of specific heats set 
to 7 = 4/3. The time-independent point-mass gravitational 
potential is $ = —GM/r, where G is Newton's constant and 
M = 1.2 Mq is the mass of the central object. The accretion 
rate ahead of the shock is M = 0.36 Mq s _1 , which is held 
fixed during the simulations. 

Our numerical simulation code, GenASiS, solves the adia- 
batic, non-relativistic, ideal MHD equations including gravity 
(cf. Eqs. (l)-(4) in Paper I). Starting from the semi-analytic 
initial condition, balance equations for mass density p, mo- 
mentum density S = pu, and magneto-fluid energy density 
Ciiuij = e to , + e kin + e mag are evolved with a second-order HLL- 
type ideal MHD scheme in a manner that ensures conservation 
of mass and energy (i.e., volume integrals of p and e„ aid + e grav ) 
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Table 1 

overview of computed models. 

B (G) loicm 2 ^ 1 ) t 


«d (ms) 


BIOLO.OEOO 
B10L1.5E15 
B10L4.0E15 


1 x 10 1U 
1 x 10 10 
1 x 10 1() 


0.0 
1.5 x 10 15 
4.0 x 10 15 


1100 

878 
678 


B12L0.0E00 a 
B12L1.5E15 
B12L4.0E15 


1 x 10 i2 
1 x 10 12 
1 x 10 12 


0.0 

1.5 x 10 15 
4.0 x 10 15 


1126 
1000 
644 


B13L0.0E00 


1 x 10 id 


0.0 


1100 



Note. — a Model computed with multiple grid resolutions. 



to numerical precision. The magnetic induction equation is 
evolved in a diver gence-free manner via the method of con- 
strained transport (Evans & Hawley 1988). (See Paper I and 
the references therein for further details. See also Appendix 
[A]in this paper.) 

Without initial perturbations the initial configuration in Fig- 
ure [T] remains stationary. In order to initiate the SASI we 
perturb the initial condition by adding small (~ 1%) ran- 
dom perturbations to the initial pressure profile in the region 
r e [-Rpnsi Rsk]- These perturbations initiate the SASI and 
allow us to study the evolution of magnetic fields in SASI- 
driven flows. 

The topology, strength and distribution of magnetic fields 
in core-collapse supernova progenitors are highly uncertain. 
A similar uncertainty applies to our knowledge of the angular 
momentum distribution in the progenitor core. These uncer- 
tainties then apply directly to the initial conditions of simula- 
tions aimed at studying the evolution and impact of magnetic 
fields in core-collapse supemovae. 

Rotation and magnetic fields in stellar interiors are in- 
timately coupled in a complex multidimensional interplay. 
Stellar core rotation can drive the evolution of magnetic fields, 
while the magnetic fields can play an important role in dis- 



tributing the core's angular momentum (e.g., Spruit 20021. 



Three-dimensional stellar evolution models (even without 
magnetic fields) extending all the way to iron core collapse 
are currently not available. However, some insight into the 
issue of core magnetic fields (and rotation) is pro vided by re- 
cent stellar evolution calculations (e.g., Heger e t al. ||2005 
Meynet et al. 201 1] >. In particular, Heger et al. (2005) in 
eluded magnetic fields in their calculations, and found that 
magnetic torques can significantly reduce the rotation rate of 
the pre-collapse iron core. The resulting magnetic fields were 
dominated by a toroidal component B v (B v /B r = 10 3 — 10 4 , 
where B r is the poloidal (radial) component of the magnetic 
field). They also reported that the core rotation rate and mag- 
netic field strength at the pre-supernova stage is an increas- 
ing function of progenitor mass. In the iron core of their 
15 Mq model, the toroidal and poloidal magnetic fields are 
B v w 5 x 10 9 G and B r w8x 10 5 G, respectively, while in 
their 35 Mq model, the toroidal and poloidal magnetic fields 
are B v sw 1 x 10 10 G and B r « 1 x 10 7 G, respectively. Ac- 
counting for the three orders of magnitude increase attained 
during core-collapse, the Heg er et al. | f2005) models predict 
the post-bounce toroidal and poloidal magnetic fields to be in 
the range of 10 12 - 10 13 G and 10 9 - 10 10 G, respectively. 
This is in the range of 'common pulsars', inferred from obser- 
vation s of pulsar spin periods and corresponding spin-down 
rates (Lorimer & Kramer 2005), but significantly lower than 
that of magnetars ( Duncan & Thompson 1992 ). 
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Investigating the role of initial _B-field topology in our sim- 
ulations is beyond the scope of this study, which is restricted 
to an initially radial (split monopole) magnetic field configu- 
ration; B r = sign(cos ■&) x Bo(R PNS /r) 2 , where d is the polar 
angle. We only vary the strength of the initial magnetic field 
B at the surface of the PNS (r = i? PNS ). In particular, we 
vary B in the range from 1 x 10 10 G to 1 x 10 13 G (cf. Table 
fTh. The initial magnetic energy density profile for the model 
with Bn = 1 x 10 10 G is represented by the dashed line in 
Figure Uj where it has been boosted by a factor of 10 6 to be- 
come visible on the plot. (The corresponding profile for the 
model with B = 1 x 10 13 G is identical to the dashed line 
in Figure [TJ) Clearly, when comparing the magnetic energy 
density to e tin and e im , all our models are initiated with weak 
magnetic fields. 

From the perspective of the Heg er et al. | ( |2005| l models, 
our initial magnetic fields are purely poloidal and stronger 
than their predicted poloidal fields, but comparable (in mag- 
nitude) to their predicted toroidal magnetic fields. Based on 
the expected multidimensional character of the post-bounce 
supernova dynamics (e.g., convection and the SASI) and the 
strength of the magnetic fields (e mag is small relative to e kin 
and e inl ) we do not expect the magnetic fields to retain the 
anisotropic (B v /B r ^> 1) configur ation p redicted by the stel- 
lar evolution calculations of Heger et al. ( 2005) 1 . Progenitors 
from multidimensional stellar evolution calculations may de- 
viate significantly from their spherically symmetric counter- 
parts (Arnett & Meakin 2011). We believe the initial mag- 
netic field configuration we have chosen has only (at best) a 
secondary impact on our results, and that initial insight into 
the MHD evolution in core-collapse supemovae can be ob- 
tained from the simulations presented here. 

For compariso n, the non-rotating simulations recently pre- 
sented by Obergaulinger & Janka | ( |20 1 1 [ > with weak initial 
magnetic fields (models sl5-B10 and sl5-Bll in that study) 
start with purely poloidal pre-collapse core magnetic fields of 
1 x 10 10 G and 1 x 10 11 G, respectively. After core-collapse 
and shock stagnation, the strength of the magnetic field in the 
stable layer separating the PNS convection zone and the gain 
region is a bout 4 x 10 12 G and 3 x 10 13 G , respectively (cf. 
Table 2 in Obergaulinger & Janka (2011 )). This layer coin- 
cides roughly with our inner boundary at r = i? PNS . Thus, the 
magnetic field strength in their collapsed weak-field models 
is comparable (initially) to that of our strongest-field model. 

Our rotating models are initiated by setting the pre-shock 
gas into rotation about the z-axis by specifying the azimuthal 
velocity according to u v = Iq sin "d/r, where Iq is the (con- 
stant) specific angular momentum. We have computed rotat- 
ing models where Iq has been set to 1.5 x 10 15 cm 2 s _1 and 
4.0 x 10 15 cm 2 s _1 . 

The discretized ideal MHD equations are solved in a cubic 
computational domain with sides L and volume Vl = L 3 . 
Cartesian coordinates are employed. The computational do- 
main is divided into N zones in each spatial dimension, re- 
sulting in N 3 cubic zones with sides Al = L/N . To conserve 
computational resources we start our simulations in a rela- 
tively small computational domain with L = L mi „ = 600 km 
and N = 512 (resulting in Al w 1.17 km). The time-step in 
our simulations (limited by the Courant-Friedrichs-Lewy con- 
dition) is about 5-10 /is, depending on the stage of the partic- 
ular run. The runs are evolved to a physical time of about 
1 s, which results in about 10 5 time-steps taken per simu- 
lation. The MHD solver is parallelized using the Message 
Passing Interface (MPI), and the computational domain is 



subdivided into blocks containing an equal number of zones, 
which are distributed among MPI processes. During the sim- 
ulations we keep the number of zones per block (MPI pro- 
cess) fixed to 32 3 . Once the SASI evolves into the nonlinear 
regime the volume encompassed by the shock grows, and the 
shock eventually interacts with the boundary of the compu- 
tational domain. When this happens, we expand the compu- 
tational domain by adding a layer of 32 3 -zones blocks (i.e., 
we add 64 zones in each coordinate direction) and restart 
the simulation from the last checkpoint written before the 
shock interacted with the boundary of the computational do- 
main. We repeat this process, and run our simulations until the 
shock interacts with the boundary of the largest computational 
box L = L mm = 1500 km, or the simulation time reaches 
t = 1100 ms, whichever occurs first. Since we keep Al fixed 
during the simulations, the largest computational domain is 
covered by 1280 zones in each spatial dimension. During a 
run, we write simulation output for analysis and visualization 
every 2 ms of physical time, resulting in tens of Terabytes of 
data from each model. 

3. SIMULATION RESULTS 

We focus on the magnetic field evolution during the non- 
linear phase of the SASI, during which magnetic fields are 
amplified most effectively and the potential for back-reaction 
of the induced fields on the fluid flow is greatest. We do 
not apply any rigorous criterion for the onset of nonlinear- 
ity; we simply note when the accretion shock deviates notice- 
ably from its spherically symmetric initial shape, and the post- 
shock velocity field has developed a significant non-radial 
component. (The upper panels of Figure [2] are representative 
of the early nonlinear phase.) 

Our simulations vary in initial magnetic field strength 
and spatial resolution, and feature both non-rotating and ro- 
tating configurations (see Table fTh. We focus first — and 
predominantly — on non-rotating models, often referring to 
model B10L0.0E00 (non-rotating model with B = 1 x 
10 10 G) as the "weak-field model," and to model B13L0.0E00 
(non-rotating model with Bq = 1 x 10 13 G) as the "strong- 
field model." The rotating models are also briefly discussed, 
but we find that the turbulent magnetic field amplification is 
mostly unaffected by rotation. 

We initiate the SASI with random pressure perturbations in 
the post-shock flow in order to avoid biased excitation of par- 
ticular modes (i.e. "sloshing" vs. "spiral"), and find that all 
our simulations exhibit flows typical of the spiral mode. This 
is consistent with ou r results in Paper I and with Blondin & 
Mezzacappa~| ( |2007] l, who found the spiral mode to dominate 
the late-time evolution independent of the ini tial pertur bation. 
It is also consistent with the conclusions of Fernan dez | ( |2010| ), 
who demonstrated that the spiral modes of the SASI can be 
viewed as a superposition of sloshing modes out of phase, 
and that any superposition of sloshing modes with non-zero 
relative phase leads to spiral modes and angular momentum 
redistribution in the post-shock flow, which potentially spins 



up the underlying PNS (see also lBlondin & Shaw |2007] >. The 
development of spiral SASI modes thus seems to be a general 
outcome (in 3D) of perturbing the (convectively stable) spher- 
ically symmetric initial condition. Moreover, a recent labora- 
tory experiment — a shallow water a nalogue to a shock insta- 
bility (SWASI, Fogli zzo et al. |2012j l — found spiral modes to 
emerge favorably from the nonlinear phase. 

3.1. Turbulence from Spiral SASI Modes 
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Figure 2. Nonlinear operation of the spiral SASI mode in model B13L0.0E00: flow Mach number |u|/cs (left panels) and magnitude of fluid vorticity 
\u>\ (= | V X u|) (right panels). The adiabatic sound speed is c$ = yyP/ p- Snapshots are taken at t = 720 ms (upper panels) and t = 820 ms (lower panels). 
To highlight the spiral mode pattern in each panel, the viewing normal is aligned with the total angular momentum in Vsh- The shock surface is traced out by the 
white contour. Velocity vectors where |u| > 10 4 km s _1 are shown in the left panels. 



Figure [2] illustrates the flows that develop from the nonlin- 
ear spiralSASI mode. The renderings are created from the 
strong-field model, but the hydrodynamic developments ex- 
hibited by this model and highlighted in Figure [2] are typi- 
cal of all our non-rotating models. The flow Mach number is 
shown in the left panels, and in the right panels the magnitude 
of fluid vorticity is displayed. 

The upper panels (t = 720 ms) depict the early develop- 
ment of the nonlinear spiral SASI mode. The shock surface 
is still quasi-spherical, but significant angular momentum re- 
distribution has occurred in the post-shock flow, and the pres- 
ence of strong counterrotating flows is app arent (cf. veloc- 



ity vectors). The shock triple-point (cf. |Blondin & Mezza- 
cappa 2007), positioned to the lower left (~ seven o'clock 



position), has just formed, and is visible as the kink in the 



shock surface (cf. white contour). The shock triple-point (a 
line segment extending across the shock surface) connects the 
pre-shock accretion flow and the two counterrotating flows 
in the post-shock gas. It moves on the shock surface in the 
counterclockwise direction in Figure [2] A layer of strongly 
sheared flows extends from the triple-point, downstream from 
the shock. This is clearly seen in both plots of fluid vorticity. 
This shear flow is one site of post-shock vorticity generation 
in our simulations. The post-shock flows are still subsonic for 
t = 720 ms. 

In the lower panels (t — 820 ms) the shock triple-point has 
completed about one and a half revolutions along the shock 
surface and is now positioned to the upper right (^two o'clock 
position). The shock volume has grown by more than a factor 
of three compared to the upper panels, and the shape of the ac- 
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cretion shock and the mass distribution in the shocked cavity 
are even more aspherical. The supersonic pre-shock accre- 
tion flow impinges on the shock at an oblique angle due to the 
aspherical shock and its off-center position. The significant 
tangential velocity component (relative to the shock surface), 
which is preserved across the shock, leads to supersonic post- 
shock flows ahead of (and directed towards) the triple-point. 
These supersonic flows, which strengthen the shear flow dis- 
cussed above, are directed down towards the PNS and result in 
further vorticity generation as they decelerate up the density 
gradient or impinge on the PNS. 

The inviscid vorticity equation is obtained by taking the curl 
of Euler's equation (e.g., Landau & Lifshitz 1959) 



du> 1 

— = V X (U X W ) + -^V/3 x VP, 

at p A 



(1) 



where the first term on the right-hand side describes changes 
in vorticity due to fluid motions, and the second term is the 
baroclinic vector. (Magnetic fields are neglected in Eq. (fTb, 
but it remains valid for weak magnetic fields.) In particu- 
lar, vorticity may be generated in regions where isosurfaces 
of density and pressure intersect. Figure [3] displays the poly- 
tropic constant k = P/p 1 (a proxy for entropy) in a slice 
through model B13L0.0E00 at t = 820 ms with focus on the 
shear flows emanating from the shock triple-point. Contours 
of constant density (dashed) and pressure (solid) are also plot- 
ted. The density and pressure contours are mostly parallel, but 
diverge strongly in the shear layer, indicating intersecting den- 
sity and pressure isosurfaces and vorticity generation through 
the baroclinic term in Eq. (fT). (Also, significant vorticity 
amplification occurs in the shear layer and elsewhere through 
the "advection" term in Eq. (fTh.) Vorticity is generated, am- 
plified, and distributed in a large fraction of the post-shock 
volume during operation of the SASI. (Movie 1 in the on- 
line material shows the generation and evolution of vorticity 
in the time interval from t = 720 ms to t = 820 ms.) The 
vorticity field exhibits strong intermittency in the late stages 
of SASI evolution. We also note the development of vorticity 
tube structures (vortex tubes) during the operation of the spi- 
ral mode. (See also Movie 2, which show s a full revo lution 
o f a v orticity still-frame at t = 820 ms.) Mee & Branden- 

tirg ( 2006 ) pointed out that the presence of vorticity may be 
:lpful for turbulent magnetic field amplification. 
Strongly forced accretion-driven turbulence develops as a 
result of the SASI, and the post-shock flow becomes roughly 
divided into a supersonic (driving) component and a subsonic 
(volume-filling) turbulent component (cf. lower left panel 
in Figure K|, This is also reflected in the probability den- 
sity function (PDF) of the velocity field below the shock. In 
the left panel in Figure |4] we plot normalized PDFs of the x- 
component of the velocity. We plot the total PDF (solid black 
line), associated with the subsonic and supersonic flow, and 
the PDF associated with the subsonic flow only (dashed black 
line). The supersonic flows contribute only to the tails of the 
distribution. The center of the distribution moves in response 
to the triple-point's movement along the shock surface (cf. 
gray curves), but when averaged over one revolution about 
the PNS, the PDF is practically centered on zero: we find 
K = Ux/«*J™ = SZo £PDF(0 dt; « 0.019. In the right 
panel of Figure HI we plot the PDF of the x-component of the 
vorticity. The vorticity PDF is more peaked, has extended ex- 
ponential tails, and is also centered about zero; [uj x /uj„ m \ FDP *=* 
0.002. Similar vorticity distributions were found in simula- 




-50 50 100 150 200 250 300 
[km] 

FigureS. Slice through model B13L0.OE0O at t = 820 ms showing the 
distribution of the polytropic constant K = Pj ' p 1 around the shear layer 
associated with the shock triple-point. Contours of constant density (p = 6 X 
10 8 g cm -3 and p = 3 X 10 8 g cm -3 ; dashed black and gray, respectively) 
and pressure (P = 1.7 X 10 27 erg cm -3 and P = 6.5 X 10 26 erg cm -3 ; 
solid black and gray, respectively) are also plotted. 



ti ons o f convectively driven turbulence by Brandenburg et al. 
(1996) and attribute d to intermittency in ny drodynamic tur- 
bulence (see also |Kraichnan || 19901 |Ishihara et al. ||2009| >. 
[Brandenburg et al. ( 1996]) characterized the intermittency of 
a variable J by the kurtosis of its PDF 



Kurt(/)=[/ 4 ] pDF /[/ 2 ] 



2 

PDF : 



(2) 



where [/"] PDF = J^ f n PDF(/) df is the rc-th moment of 
the PDF (assuming zero mean). For the x-component of the 
velocity below the shock we find Kurt('U x /w mls ) « 4.6, and 
for the ^-component of the vorticity we find K uft(u) x /u) ms ) « 
26.7 (i.e., similar to Brandenburg et al. |1996) . 

3.2. Amplification of Weak Magnetic Fields from 
Turbulence: Elementary Concepts 

The SASI-driven hydrodynamic developments result in tur- 
bulent amplification of initially weak magnetic fields, which 
is the focus of this study (see also Paper 1 1. Here we very 
briefly review some elementary concepts pertaining to such 
magnetic field amplification (for details , see for exampl e re- 
views by |Ott |1998t|Brandenburg & Subramanian |2005[ ). 

Stellar interiors are extremely good electrical conductors. 
In a perfectly conducting fluid the electric field is E = — (u x 
B), and Faraday's law (the induction equation), which gov- 
erns the evolution of the magnetic field, becomes 



dB 

~dt 



V x (u x B), 



(3) 



where the right-hand side (the induction term) describes 
changes to the magnetic field due to fluid motions. We note 
that, modulo the baroclinic vector, Eqs. ([T]) and <|3j have iden- 
tical form, suggesting a possible analogy between u> and B 



(Batchelor 1950). Moreover, Batchelor|( |1950| ) argued that 
the distribution of u> and B will be similar in fully developed 
turbulence. An important difference, however, is that the vor- 
ticity equation is nonlinear, while the induction equation is 
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Figure 4. Normalized probability density functions (PDFs) of the rc-component of velocity (left) and the ^-component of vorticity (right). The PDFs are 
constructed from the post-shock flows in model B13L0.0E00 during the nonlinear operation of the spiral SASI mode (from t = 810 ms to t = 922 ms, 
which corresponds to roughly one full revolution of the shock triple-point about the PNS). The rms values of velocity and vorticity below the shock are u rms = 
^/2E kin /M sh and w nm = y?2Q/V S h, respectively, and E kia , M sh , and C are the kinetic energy, mass, and enstrophy in V sh , the volume bounded by the shock 
surface and the surface of the PNS. We show PDFs for individual time states in gray and the average over all the time states in black. In the left panel, the total 
PDFs are represented by the solid lines, while the dashed lines are the PDFs constructed from subsonic flows only (|u|/cg < 1). The (averaged) PDF associated 
with the subsonic flow fits well with the Gaussian 0.01975 X exp [—3.1 X (u x /u rms ) 2 ]. (The PDFs constructed from the other velocity and vorticity components 
look similar.) 



linear for a specified velocity field. Nevertheless, similarities 
between vorticity and magnetic field are observed in our sim- 
ulations. 

Equation (pi can be combined with the mass conservation 
equation to form (e.g., |Landau & Lifshitz |1960| 

d /B\ _ /B 

dt W ~ v7" 



(4) 



where d/dt — d/dt + (u • V), and B/p changes due to gradi- 
ents in the velocity field. Equation |4| has an important phys- 
ical interpretation. It has the exact same form as the evolution 
equation for an infinitesimal "fluid line" connecting fluid el- 
ements and moving with the flow. Thus, two infinitely near 
fluid elements initially connected by a magnetic field line re- 
main on that magnetic field line, and the value of B/p varies 
in proportion to the distance between the fluid elements (Lan- 
dau & Lifshitz 1960). Thus, the magnetic field is "frozen 
in a perfectly conducting fluid. In an approximately incom- 
pressible fluid, the magnetic field grows in direct proportion 
to the separation between fluid elements. 

The interpretation of Eq. (HI) is equivalent to the following 
simple consideration of a magnetic flux tube, with strength b, 
length I, and cross-section a, which permeates (and is frozen 
in) a fluid element with density p: let the fluid element be 
"stretched" by the flow to a new state (characterized by b', 
V , a', and //)■ Then, mass conservation (p'l'a' = pla) and 
magnetic flux conservation (b'a' = ba) gives 



ft p 



(5) 



In the incompressible limit, the field is amplified in direct 
proportion to the stretching of the tube. At the same time, 
the flux tube undergoes a decrease in the scale perpendicu- 
lar to the stretching (a' — a x (l/V) < a). The decrease in 



flux tube cross-section proceeds until (1) the field becomes 
strong enough to react back on the fluid, preventing further 
stretching, or (2) resistive (non-ideal) effects become impor- 
tant (Section 3.7 1, or a combination of (1) and (2). Stretching 
is a very useful concept for turbulent S-field amplification. 

The frozen-in condition can result in rapid magnetic field 
amplification in a turbulent flow. Initially adjacent fluid ele- 
ments separate rapidly, perhaps exponentially, in the chaotic 
flows that characterize turbulence (Ott 1998). Thus, an ini- 
tially weak magnetic field (i.e., u is independent of B) may 
amplify exponentially by stretching, and the growth rate is 
roughly given b y the inverse t urnover time of turbulent ed- 
dies (e.g., Kulsrud & Anderso n |1992[ ). This is also apparent 
from a simplistic consideration of Eq! (HJi: the turbulent ve- 
locity varies significantly ~ 0(u'™J over a turbulent eddy of 
size A eddy , hence B~ 1 (dB/dt) ~ u'™ /A rfdy . Exponential am- 
plification of weak magnetic fields is commonly seen in MHD 
turbulence simulations (e.g., |Cho & Vishniac 2000 1 Branden- 
burg 2001; Haugen et al. 2004). Exponential growth ceases 
when the magnetic field becomes strong enough to cause a 
back-reaction on the fluid (i.e., u becomes dependent on B). 

Amplification of weak magnetic fields through turbulent 
stretching is initially a kinematic mechanism (i.e., described 
by Eq. p) for a specified velocity field). As such, it differs 
from magnetic field amplification by the MRI in a fundamen- 
tal way. The MRI is a dynamic mechanism, described by the 
full MHD system of equations, and requires the Lorentz force 
to be included in Euler's equation. (Also, the MRI requires 
differential rotation in the PNS to operate, while a turbulent 
dynamo can operate without PNS rotation.) However, for 
weak progenitor i?-fields (Heger et al. 2005[ ), both mech- 
anisms req uire high spatial resolut ion for simulation (e.g., 
Obergaulinger et al. 2009] |Paper I| l, and may ultimately be 
computationally prohibitive to capture properly in large-scale 
multi-physics simulations. 
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We comment here on the amount of kinetic helicity u> ■ u 
in our simulations of SASI-driven turbulence. Kinetic he- 
licity is a measure of "handedness" (or lack of mirror sym- 
metry) in the turbulent flows, and is an important quantity 
in dynamo theory (e.g., Brandenburg & Subramanian 2005, 
and reference therein). Turbulent flows with kinetic helic- 
ity can support a so-called inverse cascade and produce large 
scale magnetic fields (i.e., larger than the turbu lent forcing 
scale) in the nonlinear, s aturated state (e.g., Meneguzzi et 
[aT~|[T98T] |Brandenburg~||200T] l . Non-helical turbule nce re- 
sults in mo stly sma ll-scale magnetic fields (e.g., Branden- 
burg et al ~| [1996 ; Hau gen et al. ||2004] >. We have con- 
structed PDPs of the relative kinetic helicity h km = (u) ■ 
u)/(w mls w, ml! ) in the post-shock flow. The kinetic helicity dis- 
tributions are similar to the vorticity distributions (strongly 
peaked with exponential tails), with [/i ki „] PDF ~ —8.3 x 10~ 4 
and Kurt(/i kin ) ss 28.4. Despite the apparent handedness as- 
sociated with the spiral SASI mode, the resulting turbulence 
is essentially non-helical. (This may change, however, if a 
rapidly — and differentially — rotating PNS is included in the 
model.) Thus, we expect that SASI-driven turbulence in our 
simulations results in magnetic field amplification due to a 
non-helical small-scale dynamo. 

3.3. Time Evolution of Global Quantities 

An overview of the simulations with non-rotating initial 
conditions is given in Figure BJ in which we plot the time- 
evolution of selected globally integrated quantities for models 
B10L0.0E00 (solid lines), B12L0.0E00 (dashed lines), and 
B13L0.0E00 (dotted lines): kinetic energy (upper left), rela- 
tive magnetic energy change (upper right), angular momen- 
tum (lower left), and average shock radius (lower right). All 
quantities are obtained from integration over the shock vol- 
ume V sh , bounded by the surface of the PNS dV ms and the 
surface of the shock dV sh . The magnetic energy change is 
scaled with the initial magnetic energy for easy comparison 
across the models ((E mK — E m . df , ) /E^ B is plotted). 

The kinetic energy of the settling flow beneath the shock is 
initially about 2 x 10 -3 B. It begins to grow rapidly during the 
initial ramp-up phase of the SASI, which starts around 400 ms 
in all models. In particular, for the weak-field model, the 
post-shock kinetic energy grows exponentially with a nearly 
constant growth rate over the time period extending from 
t « 510 ms to t sa 720 ms. The growth time during this 
epoch is about r ~ 85 ms (cf. dash-dotted line in the top 
left panel in Figure BJ). The kinetic energy in the models 
with a stronger initiafmagnetic field grows somewhat slower 
initially (t < 660 ms), and then at a rate similar to that of 
the weak-field model. The growth slows down considerably 
for t > 800 ms, but the kinetic energy continues to grow 
throughout the nonlinear phase and reaches similar levels in 
all three models, with variability on a shorter timescale su- 
perimposed. When averaged over the time interval extend- 
ing from 900 ms to 1100 ms we fincf^-EL)^ 8 = 0.051 B, 
0.048 B, and 0.044 B for models B10E0.0E00', B12L0.0E00, 
and B13L0.0E00, respectively. While these time-averaged 
post-shock kinetic energies are slightly smaller in the mod- 
els with a stronger initial magnetic field, we have not found 
convincing evidence that this slight decrease in global kinetic 
energy is a result of a stronger magnetic field. We find, how- 

7 The temporal average of a variable X, over the interval T = fej — tl, is 
denoted {X)\\ = A f* 2 X dt. 



ever, that strong magn etic fields affect flows on small spatial 
scales (cf. Section BCo} . 

In terms of spherical harmonics, the SASI is charact erized 
by exponentially grow ing power in low-order modes ( |Blondin| 
& Mezzacappa 2006). As a result of this, the shock sur- 
face deviates exponentially from its initially spherical shape. 
The obliquity of the shock front relative to the pre-shock ac- 
cretion flow causes the non-radial post-shock kinetic energy 
to gr ow (also exponential ly) at the expense of thermal en- 
ergy (Blon dmet al.|[2003] l. We have decomposed the post- 
shock kinetic energy into radial and non-radial components; 
^;„,|| = I J Vsh pu 2 r dV and E klIh± = \ J^ p{u% + u%) dV, 
respectively. The non-radial component grows much faster 
(t < 50 ms) than the radial component, and the growth seen 
in Figure BJ is due to a combination of the two. (The kinetic 
energy associated with the three velocity components become 
similar in the saturated state, and _E kin .j_ ~ 2_E kin || .) Saturation 
of the post-shock kinetic energy may be due to the develop 



mentof turbulence via secondary instabilities (e.g., Guilet et 



al. 2010 and Section 3.6 in this paper). The exact details that 



determine the growth rate of the post-shock kinetic energy are 
tied directly to the physical origin of the SASI, which is not 
the focus of this paper. We focus primarily on magnetic field 
amplification in the flows that result from SASI activity. 

The magnetic energy grows at the expense of th e turbulent 
kinetic energy below the shock (cf. Sectio n |3.6| l. After an 
initial spurt, all the models shown in Figure pi experience an 
early period of exponential magnetic energy growth with es- 
sentially the same growth rate (cf. the temporal window from 
t = 650 ms to t — 780 ms). Such evolution is expected in a 
kinematic growth regime, in which the magnetic field's back- 
reaction on the fluid is negligible. The magnetic energy in the 
weak-field model (B10L0.0E00, B a = 1 x 10 10 G) grows 
exponentially at a nearly constant rate, with growth time 
r sa 66 ms, until the end of the simulation (t = 1100 ms), and 
receives a total boost of about five orders of magnitude. The 
magnetic energy growth time is significantly shorter (~ 25%) 
than the total kinetic energy growth time during the over- 
lapping period of exponential growth. In the model with 
B = 1 x 10 12 G (B12L0.0E00) we also find that E m . d , grows 
steadily until the end of the run (t — 1126 ms). The magnetic 
energy in this model initially grows at the same rate as the 
weak-field model, but its growth rate clearly tapers off at later 
times (t > 900 ms). The strong-field model (B13L0.0E00, 
Bq = 1 x 10 13 G) also exhibits exponential magnetic energy 
growth (r s=s 66 ms) early on. Then, around t w 780 ms, 
its growth rate drops almost discontinuously, and E, m[! , grows 
by only about 50% for the remainder of the simulation (until 
t = 1100 ms). Model B13L0.0E00 receives a total boost in 
magnetic energy of about a factor of 300. The abrupt change 
in the magnetic energy growth rate observed in the strong- 
field model occurs when magnetic fields become dynamically 
imp ortant in localized regions below the shock (cf. Section 

At the end of the simulations the magnetic energy in mod- 
els B10L0.0E00, B12L0.0E00, and B13L0.0E00 has reached 
about 1.8 x 10~ 7 B, 2.3 x 10" 4 B, and 8.9 x 10~ 4 B, respec- 
tively. The magnetic energy in the weak-field model is many 
(~ five) orders of magnitude below the post-shock kinetic en- 
ergy at this point. In the strong-field model it saturates below 
10 -3 B, which is also significantly less than the total kinetic 
energy in the post-shock flow. Also, the boost in magnetic en- 
ergy in the strong-field model is almost an order of magnitude 
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Figure 5. Time-evolution of global quantities integrated over the shock volume V$ h , bounded by the surface of the PNS, <9Vpns, and the surface of the shock, 
c9Vsh, in non-rotating models in which the initial magnetic field has been varied from 1 X 10 10 G to 1 X 10 13 G. Plotted are kinetic energy E kia (top left), magnetic 
energy change (relative to the initial) AE mlg /E m ^.o (top right), angular momentum |L| (bottom left), and average shock radius R S h = (3Vsh/(47r)) 1 ' 3 . Models 
B10L0.0E00, B12L0.0E00, and B13L0.0E00 are represented by solid, dashed, and dotted lines, respectively. The initial magnetic energy -E mag ,o in these models 
is 2.3 X 10~ 12 B, 2.3 X 10 — 8 B, and 2.3 X 10~ 6 B, respectively. The dash-dotted lines in the top panels are proportional to exp (t/r), where the growth times 
are t = 85 ms and r = 66 ms in the top left and top right panels, respectively. 



lower than the difference in average post-shock kinetic energy 
between the strong-field and the weak-field models, which is 
about 7 x 1CP 3 B. At the end of the simulations the magnetic 
energy in the strong-field model is less than a factor of four 
larger than in model B12L0.0E00, while it initially was a fac- 
tor 10 2 larger. We point out that the magnetic energies listed 
above are values recorded when the simulations were stopped 
after roughly one explosion time (t ~ 1 s). (SASI-induced 
magnetic field amplification ceases once the explosion is ini- 
tiated.) However, the listed values should not be interpreted 
as upper limits on the magnetic energy for the different initial 
magnetic fields. The magnetic energy growth rate is (for rea- 
sons we detail in this paper) underestimated by the numerical 
simulations, and we suspect that all models — independent of 
the initial magnetic field — will reach a saturated state, similar 
to the strong-field model, well within an explosion time. The 



issue of magnetic energy growth, saturation, and its effect on 
the post-shock flow will be discussed in later sections. 

The induced magnetic fields do not notably affect the global 
characteristics of the shock evolution. The plots of total an- 
gular momentum |L| in V^ h and the average shock radius 
.ff sh show that these quantities reach similar values in all the 
non-rotating models. The angular momentum reaches a few 
xlO 47 g cm 2 s _1 , consistent with the comparable models of 



Blondin & Mezzacappa ( 2007) and Fernandez (2010). More- 
over, during the nonlinear evolution, after the period of ex- 
ponential growth of the angular momentum in V sh , we find 
|L| < fMRi, with / w 0.25 (cf. |Fernandez |2010| . The av- 
erage shock radius exhibits significant variability, and briefly 
exceeds 500 km in some of the models. In particular, we 
find (Rs h )o.ll = 484 km, 466 km, and 438 km, for models 
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B10L0.0E00, B12L0.0E00, and B13L0.0E00, respectively. 

The larger-amplitude oscillations in kinetic energy and 
the smaller average shock radius exhibited by the strong- 
field model (B13L0.0E00) may be attributed to this model's 
somewhat different nonlinear evolution, which is due to the 
stochastic nature of the nonlinear SASI rather than to the 
stronger initial field. Model B13L0.0E00 first evolves into 
a typical spiral mode pattern (cf. Figure [2|, but later devel- 
ops a flow pattern reminiscent of the sloshing mode, with two 
oppositely directed high-speed streams, emanating from op- 
posite sides of the shock, terminating on opposite sides of the 
PNS, or colliding head-on deep in the shocked cavity. The ap- 
pearance of this flow pattern coincides with the turnover in the 
angular momentum seen in the lower left panel (t sa 950 ms). 
Note that there is little or no response in the magnetic energy 
evolution due to these rearrangements in the large scale flow. 
This is consistent with the magnetic field being amplified by 
small-scale rather than large-scale flows. 

Comparing Figure B] with Figure 12 of Paper I we note 
that spatial resolution (A/ = 1.17 km in this paper versus 



A^ = 1.56 km in Paper I) affects global magnetic quanti 
ties much more than global fluid quantities. In particular, 
models 3DB10Rm and 3DB12Rm in Paper I correspond to 
models B10L0.0E00 and B12L0.0E00 respectively. The in- 
creased spatial resolution has no significant impact on the 
post-shock kinetic energy, total angular momentum, average 
shock radius, or exponential growth time for the kinetic en- 
ergy (r « 85 ms). The magnetic energy in model 3DB12Rm 
is boosted by a factor of 2 x 10 3 , while model 3DB10Rm 
received a boost of less than a factor of 10 3 . We also mea- 
sured an exponential growth time for the magnetic energy 
of r « 71 ms over several hundred milliseconds in Paper 
I These results are somewhat different than the results pre- 
sented in Figure [5] which show that the increased resolution 
results in a larger boost in the magnetic energy and shorter 
exponential growth time (r sa 66 ms), and that the magnetic 
energy growth at late times depends on the initial magnetic 
field strength. The sensitivity to spatial resolution was also 
pointed out in our previous study, and it will be further dis- 
cussed later in this paper. 

3.4. Evolution of Spherically Averaged Radial Profiles and 
Saturation of Magnetic Energy 

The magnetic energy in the strong-field model reaches satu- 
ration relatively early in the nonlinear evolution (t ss 780 ms). 
To help elucidate the physical conditions under which the 
magnetic energy growth in our simulations is quenched, we 
plot, in Figure[6] spherically averaged radial profiles from the 
evolution of the weak-field model (B10L0.0E00, left panels) 
and the strong-field model (B13L0.0E00, right panels). In 
the upper panels we plot the rms magnetic field strength B, ms 
(solid lines) and the maximum magnetic field strength B ma 
(dashed red lines) versus radial distance from the center of the 
PNS. Values are computed in spherical shells centered on the 
origin of the computational domain r = \J x 2 + y 2 + z 2 — 0, 
with thickness SL = 20 km and volume 5Vi, which includes 
all zones with radial coordinate r € [n — 6L/2, r*j + SL/2), 
with n = 50 km, 70 km,. . .,(L — 8L)/2. The rms magnetic 
field is computed from the shell- volume-averagecP] magnetic 
energy density B m = y/2(ia(e m ,}sv i , while B ma is simply 

8 (X) v = -h J v X dV denotes the volume average of X over the vol- 
ume V. 



the maximum magnetic field over all zones in each shell. In 
the lower panels we plot the shell-volume-averaged magnetic 
energy density (e mit .)sVi (solid lines), kinetic energy density 
(ekin)i5Vi (dotted blue lines), and fluid pressure (P)sVi (dashed 
red lines) versus radial distance. For each model we plot four 
profiles (time states) of each variable. The time states, which 
are also indicated with diamonds (B10L0.0E00) and squares 
(B13L0.0E00) in the upper right panel of FigureBJ are chosen 
to emphasize the temporal magnetic field evolution in each of 
the models and to contrast the two models. 

The spherically averaged radial profiles further illustrate 
the differences in magnetic field evolution of the weak- 
field and strong-field models. The magnetic field in model 
B10L0.0E00 intensifies steadily throughout the nonlinear 
evolution, at all radii below the shock, and B nm evolves self- 
similarly during the later stages. Toward the end of the weak- 
field run, the rms magnetic field has received a boost of about 
two orders of magnitude near i? PNS , and stays above B for 
r < 500 km. The maximum magnetic field is typically an 
order of magnitude above the rms magnetic field, which is 
an indication of strong spatial intermittency in the magnetic 
field (see also Figures|7]and[8]below). In the lower left panel 
of Figure [6] we see that the magnetic energy density is many 
orders of magnitude smaller than the kinetic energy density 
and pressure for the shown times, consistent with kinematic 
magnetic field growth. Not even in localized regions does the 
magnetic energy density become comparable to the kinetic 
energy density or the pressure. At the end of the simulation, 
there are only a few zones in which the ratio of magnetic-to- 
kinetic energy exceeds 10~ 2 . 

The strong-field model's magnetic energy evolution is not 
governed by kinematic growth (except for during the initial 
boost received at early times), but rather by dynamic interac- 
tions with the fluid in a sa tura ted s tate (and numerical dis- 
sipation; see Sections 3.6 and 3.7 1. The magnetic energy 
in this model falls off the exponential growth curve around 
t = 780 ms. The thin solid curve in the upper right panel in 
Figure|6](t = 720 ms) represents the transition from the initial 
magnetic field profile (dash-dot curve) to the saturated state. 
Although the post-shock flow is governed by vigorous turbu- 
lence at later times, there are only minor changes to the rms 
and maximum magnetic field profiles. The relative boost of 
B„ and B max in model B13L0.0E00 is about an order of mag- 
nitude less than what is observed in model B10L0.0E00. At 
the end of the strong-field model, the rms magnetic field ex- 
ceeds 10 13 G for r < 225 km, while the maximum magnetic 
field exceeds 10 14 G out to r « 200 km. 

At the end of both runs the rms magnetic field follows 
an exponential decrease with radius, B m , is proportional to 
exp (— t/Lb), where the characteristic length scale Lb is 
about 120 km and 100 km for model B10L0.0E00 and model 
B13L0.0E00, respectively (cf. dash-dotted lines). From the 
upper panels in Figure |6j it is apparent that the exponential 
decrease in B nm with radius holds reasonably well in both 
models throughout the nonlinear evolution. Moreover, the av- 
eraged kinetic energy density below the shock has increased 
significantly compared to the initial condition (Figure fTJ and 
roughly follows a power-law in radius r a , with the exponent a 
varying between —2.7 and —2.3. The decrease in kinetic en- 
ergy density with radius is mostly due to the decrease in mass 
density: the mass density falls off as r~ 3 inside r — 150 km, 



and almost as r ~ outside 



150 km. During the runs, 



the spherically averaged pressure remains relatively quiescent 
inside r = 150 km, where it falls off as r -4 . 
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Figure 6. Spherically averaged radial profiles at selected times for the non-rotating weak-field (left panels) and the strong-field (right panel) models. Upper 
panels: rms magnetic field strength (solid black) and maximum magnetic field strength (dashed red). For reference, the initial magnetic field profile is plotted 
in each panel (dash dot). We also plot reference lines proportional to cxp (—r/Lg) (dash-dotted) with Lg = 1 20 km (weak-field case) and Lg = 100 km 
(strong-field case). The magnetic field strengths have been normalized to the initial value Bn at r = Rpus- Lower panels: magnetic energy density (solid black), 
kinetic energy density (dotted blue), and thermal pressure (dashed red). For model B10L0.0E00 profiles are plotted for 666 ms, 792 ms, 918 ms, and 1100 ms 
(marked by diamonds in the upper right panel of Figure[5j, while for model B13L0.0E00 we plot radial profiles for times 720 ms, 820 ms, 922 ms, and 1100 ms 
(marked by squares in the upper right panel of Figure pi. Thicker lines indicate a more advanced time state. Note that the magnetic energy densities in the 
weak-field model (lower left panel) has been multiplied by a factor of 10 4 . 



As noted above, the decrease in kinetic energy density (the 
source of magnetic energy) follows a power law with radius, 
while the magnetic energy density (and B ms ) decreases ex- 
ponentially with radius. On the other hand, we find that the 
enstrophy ^(u! 2 )sVi also decreases exponentially with radius, 
with a length scale comparable to (but somewhat shorter than) 
that of (e-„- W. . Mo reover, by comparing the lower p anels of 
Figure 6 in Paper I (showing |B|) with Figure 11 in |Paper 1 



(showing |<_>|) we see that the spatial distribution of magnetic 
field and vorticity is very similar. These observations support 
a similarity between vorticity and magnetic field (Batchelor 
1950). (We have not investigated the physical reasons for the 
particular spatial distribution of vorticity and magnetic field 
in further detail, but we plan to do so in a future study.) 



The relative boost in _B rms decreases monotonically with the 
initial field strength in our models. The results from model 
B12L0.0E00 (not shown) confirm this trend. This is because 
the models with stronger initial fields reach saturation during 
the simulation. (The growth rate is the same in all models dur- 
ing the kinematic regime.) Saturation occurs when the mag- 
netic energy becomes comparable (locally) to the kinetic en- 
ergy. In particular, we find that the kinematic growth regime 
ends when /5 , ^ lx /(2/_ ) < (e kia )sVi- 

For model B12L0.0E00 we find BL/( 2 Mo) « 0.1 x 
(e m ) SVi and B^J(2fj, ) w 0.3 x (e kin ) SVz for t = 966 ms 
and t = 1126 ms, respectively, which represent time states 
after the magnetic energy growth has fallen off the exponen- 
tial curve with growth time t = 66 ms (Figure E\. In model 
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B13L0.0E00, the ratio B^ x /(2fj, {e an )sv i ) stays above 0.3 for 
the three most advanced time states shown in Figure |6| and 
hovers around unity at the end of the simulation. In both 
models the ratio stays remarkably constant with distance from 
the PNS, varying by less than a factor of two, out beyond 
r = 300 km (although exact details vary by model). Thus, 
turbulence-induced magnetic fields may impact dynamics in 
localized regions throughout the shock volume. 

3.5. Full Spatial Distributions, Intermittency, and 
Saturation of Magnetic Energy 

Figure 17] shows volume renderings of the magnetic field 
magnitude in model B10L0.0E00 near the end of the simula- 
tion (t = 1068 ms). The left panel shows a global view of the 
amplified magnetic fields below the shock (white contour), 
and the right panel shows a zoomed in view in a 200 km 3 
volume near the PNS. These renderings illustrate the compli- 
cated, highly intermittent magnetic fields that develop from 
SASI-induced turbulent flows. The magnetic energy is con- 
centrated in thin, folded, intense flux ropes. Notice in the left 
panel that amplified magnetic fields do not extend all the way 
to the shock, but are confined to a smaller volume, which is 
characterized by highly turbulent flows (see also the distribu- 
tion of vorticity for model B13L0.0E00 in the volume render- 
ings in Figure [2] and the associated movies). The spatial dis- 
tribution of the magnetic field (in particular the intermittency) 
in the strong-field model is similar to that of the weak-field 
model. 

Normalized PDFs of individual components of the mag- 
netic field from model B10L0.0E00 are plotted in the left 
panel of Figure [8] The shape of the distributions is strongly 
peaked with extended exponential tails (similar to the vorticity 
distributions in Figure 0), and the PDFs of the different mag- 
netic field components are practically indistinguishable. The 
intermittency is high, Kun(B x /B, m ) « Kun(B y /B mi ) « 
Kurt(B z /B nm ) ~ 32.5 (somewhat larger than, but comparable 
to, the intermittency of the vorticity field), which is consistent 
with the visual impression given in Figure |7j We note that 
the PDFs are highly symmetric, which implies small over- 
all polarity of the field; the magnitude of the mean values 
[Bv/BJ^ [By/BJl^ and [B,/BJ[„ are all less than 
0.01. The PDF of the cosine of the angle between oj and B 
is plotted in the right panel of Figure [8] which shows that the 
vorticity and magnetic field tend to be aligned or anti-aligned, 
and gives further support to the similarity between the mag- 
netic and vorticity fields. Similar distributions were also re- 
ported by Brandenburg et al. ( 1996| ). 

The ratio B^. a / {'Zji^e^) svi) (used in Section 3.4 to char- 
acterize saturation of the magnetic energy) is still only an ap- 
proximate measure of the relative strength of the magnetic 
field since the kinetic energy density is averaged over the 
shell. The highly intermittent magnetic fields cre ated by tur- 
bulence have been "expelled" from the fluid (Thompson & 
Duncan 1993| l. As evolution proceeds, an increasing percent- 
age of the total magnetic energy resides in regions where the 
ratio of magnetic-to-kinetic energy /3 k r 1 (= u 2 /|u| 2 ) exceeds 
10~ 2 , 10 _1 , and 1: 55%, 10%, and 0.5%, respectively, for 
model B12L0.0E00 at t = 966 ms, while for t = 1126 the re- 
spective percentages have increased to 72%, 20%, and 1.5%. 
(The percentages quoted for t = 1126 ms are very similar to 
those quoted in Paper I for model 3DB12Ah, which was com- 
puted with the same resolution and initial condition as model 
B12L0.0E00, but with a different initial perturbation.) The 
fraction of the total magnetic energy concentrated in regions 



with high f3 kiT } stays relatively constant during the saturated 
state of the strong-field model (t > 780 ms): we find ~ 90% 
(A" 1 > 10- 2 ), ~ 50% (ft- 1 > 10- 1 ), and ~ 7% (j3~ l > 1). 

Volume renderings in Figure|9]show the spatial distribution 
of /3" 1 late in the evolution orthe strong-field model, illus- 
trating the spatial and temporal intermittency of turbulence- 
induced strong magnetic fields. The snapshots are tempo- 
rally separated by 10 ms , wh ich is longer than the turbulent 
turnover time (cf. Section [3T6| ), but is significantly shorter than 
both the advection time and the Alfven crossing time (defined 
loosely as i? sh /|u| and R sh /v^ respectively). Concentrations 
of high f3~J-, which can briefly exceed unity in localized re- 
gions, are scattered throughout the shock volume. As noted 
above, about 50% of the total magnetic energy resides in re- 
gions where /3" 1 > 10 _1 , and these magnetic fields occupy 
less than 10% of the total shock volume. (Movie 3 and Movie 
4 in the online material show the time-evolution of ft" 1 from 
t = 1050 ms to t — 1100 ms and a full revolution of ft" 1 for 
t = 1100 ms, respectively.) 

Alfven waves are likel y excited by the SASI activity dis- 
cussed in this paper. Suzu ki et al. | d2008) performed 
simulations in spherical symmetry and investigated the role 
of Alfven waves on the core-collapse supernova explosion 
mechanism. They found that — for sufficiently strong mag- 
netic fields (> 2 x 10 15 G) — heating associated with Alfven 
wave energy dissipation may revive the stalled shock. For 
weaker magnetic fields no shock revival was observed. We 
have not attempted to identify Alfven waves, or energy dissi- 
pation due to Alfven waves, in our simulations (the highly dy- 
namic nature of the S ASI-driven flows makes this a nontrivial 
task). However, the magnetic fields attained in our strong- 
field model are significantly weaker than 10 15 G and we do 
no t expect Alfven wav e heating due to the mechanism studied 
by Suzuki et al. | ( |2008] l to result in significant energ y deposi- 
tion near th e shock (or elsewhere) in our simulations. Guilet et| 
al. ( 201 1] > recently suggested a mechanism for magnetic field 
amplification in the vicinity of an Alfven surface (i.e., where 
v A = |u|). In their model, Alfven waves, excited for example 
by the SASI, may amplify near the Alfven surface and create 
a dynamic back-reaction. However, the turbulent nature of the 
hydromagnetic evolution in our simulations may result in un- 
favorable conditions for this mechanism to operate. In partic- 
ular, regions where v A — |u| appear and disappear in a highly 
intermittent manner (cf. Figure 9\ and associated movies). 

3.6. Spectral Analysis 

Further important insight into the numerical simulations 
can be gained from a Fourier decomposition of the magnetic 
and k inetic energy. In particular, our analysis presented in Pa- 
per I lacked the ability to quantify the amount of turbulent ki- 
netic energy available to amplify the magnetic field as well as 
the magnetic field's impact on the evolution of the small-scale 
flows. We se ek to address these questions in this section. 

Following Ryu et al. (2000) we compute Fourier ampli- 
tudes from components of the velocity and magnetic fields in 



the computational domairF] 



X(k) = 1 / X(x) 

v L J Vl 



exp (zk • x) dV, 



(6) 



where X(x) represents ^/puj or Bj, with j £ {x, y, z}. We 

9 The Fourier transforms are computed using the FFTW library docu- 
mented at http://www.fftw.org 
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Figure 7. Magnetic field magnitude |B| near the end (t = 1068 ms) of the weak-field simulation (B10L0.0E00). The left panel shows a global view of the 
magnetic fields below the shock (traced out with a white contour). The right panel shows the magnetic field in a 200 km 3 volume near the PNS. 
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then compute the kinetic and magnetic spectral energy density 
on a /c-space shell 

£»(*) = £/ V|Vwl 2fc2 ^fe ( 7 ) 

1 Jk-shell 



and 



,(*) 



1 
2^o 



fc-shell 



3 



BA 2 k 2 dQ k , 



(8) 



respectively. The magnitude of the wave vector (wavenum- 
ber) is k = |k| = (k 2 , + k 2 y + k 2 z ) 1 / 2 and dQ k is a solid 
angle element in Fourier space. Proper normalization of the 
Fourier transform ensures that integration of the spectral en- 
ergy densities over fc-space equals real-space integrals of the 



corresponding energy densities; i.e., 



and 



e ti „ dk= e kin dV = E kl 
fc„,i„ Jv L 



e milg dk = I e, m dV = E m 



Vl 



(9) 



(10) 



where the integrals over k extend from fc mi „ = 2ir/L (de- 
fined by the spatial scale of the computational box L) to 
fc max = 2ir/Al (limited by finite grid resolution). The real- 
space integrals extend over the volume of the computational 
domain Vl- Note that _E ma „, the total magnetic energy in the 
computational domain, practically equals the magnetic energy 
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below the shock wave E m ., s , while a similar equality does not 

hold for the kinetic energy (E ti „ > E m ) because the kinetic 
energy of the supersonic accretion flow above the shock (in- 
cluded in the Fourier transform) is substantial, but contributes 
mostly to the spectrum at small fc (see Figure 12 1. 



3.6. 1 . Varying the Initial Magnetic Field Strength 

The evolution of magnetic energy in Fourier space is shown 
in Figure [lOJ in which we plot magnetic energy spectra dur- 
ing the nonlinear SASI for the weak-field model (left panel) 
and the strong-field model (right panel), for the same times 
used for the spherically averaged radial profiles displayed in 
Figure[6] Initially (not shown), the magnetic energy spectrum 
decreases monotonically with increasing fc. Then, most of the 
magnetic energy resides on relatively large scales (k < 0.1), 
and e mag is roughly proportional to fc~ 2 for larger fc-values. 
The energy spectra in Figure 10 exhibit features typical of 



MHD turbulence simulations (see Brandenbu rg & Subrama- 
nian |2005, for a recent comprehensive review). The spectral 



magnetic energy density increases with wavenumber, roughly 



as e„, 



fc 3 / 2 , for small fc-va lues (cf. JKazantsev 



1968 



Brandenburg & Subramanian 2005). It reaches a maximum 
around 0.2-0.3 km -1 , where the turnover is due to numeri- 



cal diffusivity (Figure 14 1, beyond which it decreases rapidly 
with increasing k. 

The spectral magnetic energy density remains basically 
self-similar over the time-intervals displayed in Figure [10] 
The weak-field model exhibits the exponential growth on all 
scales typical of a kinematic small-scale dynamo (Branden- 
burg & Subramanian |2005) , the peak value increasing by al- 
most three orders of magnitude (from ~6x 10~ 10 B km to 
~4x 10~ 7 B km), as seen in Figure [lO] In contrast the peak 
value saturates in the strong-field modek increasing by only a 
factor of about four (to ^ 2 x 10~ 3 B km). However the spec- 
tral shape stays relatively unchanged in both cases, with the 
full width at half maximum being roughly constant over time 
(~ 0.40 km -1 and ~ 0.38 km -1 for the weak-field model 
and the strong-field model, respectively). 

The spectral shape in Figure [10] appears unchanged across 



the two models, even though the strong-field model saturates 
and the weak-field model does not. There are, however, small 
differences. At the end of the simulations (t = 1100 ms), the 
normalized spectra ^ mag /max(e mag ) of models B10L0.0E00 
and B12L0.0E00 lie practically on top of each other for all 
k. The corresponing spectrum of the strong-field model fol- 
lows those of the weaker-field models for large fc (although 
the peak is slightly shifted to the left), but has "excess" power 

f or fc < 0.1 km -1 : the integral L — r f, ' e m ,„, dk is about 

^ fe max(e mag ) J k m [ n m ^ 

65% larger in the strong-field model than in the weaker-field 
m odels. In the simulations of non-helical MHD turbulence 
by |Hau gen et al. (2004), the magnetic energy spectra grow 
self-similarly until saturation (e mag ~ e kin ), which occurs first 
on smaller spatial scales, and on larger scales later (i.e., the 
smallest wavenumber where e mag ~ e a „ moves to even smaller 
fc), and the magnetic energy spectrum appears to align itself 
with the kinetic energy spectrum with e mag > e ki „, almost up 
to the forcing scale. This may suggest that the shape of the 
saturated and unsaturated spectra should differ more than dis- 
played in Figure 10 Despite the differences (saturated or not), 



we have not been able to find a reasonable explanation for 
why the spectral shape in the two models remain so similar. 

From the spectral magnetic energy distribution we obtain 
the mean magnetic wavenumber 



1 
E m 



fc6 mag dk, 



(11) 



and the characteristic spatial scale of the magnetic field A nm „ = 
27r/fc mag . (The mean magnetic wavenumber is indicated by a 
diamond on each of the energy spectra in Figure 10 ) During 



the initial ramp-up to nonlinear SASI evolution we find that 
A mag decreases rapidly with time, from A m , le s=s 60 km initially 
to A m , le w 20 km around t = 650 ms, and stays relatively 
constant thereafter in all the non-rotating models. Magnetic 
field amplification in our simulations is caused by turbulent 
stretching of flux tubes. In a kinematic dynamo the charac- 
teristic scale of the magnetic field decreases exponentially. 
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Figure 10. Temporal evolution of the spectral magnetic energy density e mli „ for the weak-field (left) and the strong-field (right) models, respectively. The spectral 
distributions are plotted at 666 ms, 792 ms, 918 ms, and 1100 ms (model B10L0.0E00) and 720 ms, 820 ms, 922 ms, and 1100 ms (model Bf3L0.0E00), 
respectively (i.e., the same times for the respective models as those displayed in Figure[6). The dotted vertical reference lines denote spatial scales (from left to 
right) of 300 km, 20 X Al, and 10 X Al, w here Al = 1.17 km is the size of a computational cell. The dash-dot line in each panel is proportional to fc 3 ' 2 . The 
mean magnetic wavenumber fc mag (Eq. Ill)) is indicated by diamonds. 



If the kinematic approximation remains valid, the decrease 
is halted by resistive dissipation when the spatial dimension 
of the field (the flux tube thickness) approaches the resistive 
scale (Schekoch ihin et al. |2002 1. The kinematic approxima- 
tion remains valid throughout the evolution of the weak-field 
model. 

The temporal constancy of A rallg for t > 650 ms in the weak- 
field model is a strong indication that numerical diffusion 
plays an important role in our simulations. In our nu merical 
scheme, we have adopted the HLL Riemann solver ( |Harten| 
et al. |1983]l, which approximates the MHD Riemann prob- 



lem by only considering the left- and right-propagating fast 
magnetosonic waves. This approximation results in diffusive 
evolution of intermediate waves (e.g., slow magnetosonic, 
Alfven, and entropy waves), and is the main source of dissipa- 
tion in our simulations. (No other form of dissipation, phys- 
ical or numerical, has been explicitly included in our simula- 
tions.) The inherent diffusivity of schemes based on the HLL 
Riemann solver also affects the evolution of small-scale struc- 
tures (e.g., turbulence induced magnetic fields. See Appendix 
lAl for further details on the source and nature of numerical 
dissipation of the magnetic energy in our simulations.) 

Moreover, in the strongly nonlinear regime of the SASI at 
t > 750 ms we find that A mag is somewhat larger (~ 10%) 
in models with a stronger initial magnetic field. Specifically, 
we find (A mag )Jg' « 18 km for the weak-field model and 
(•\» ag )o'8s ~ 20 km for the strong-field model. This trend is 
consistent with the magnetic field becoming strong enough to 
cause a back-reaction on the fluid through the magnetic ten- 
sion force and thereby limit the extent to which magnetic flux 
tubes are stretched and bent by the chaotic flow induced by 
the SASI. Schekochi hin et al~| ( 2001| l observed a strong anti- 
correlation between the magnetic field strength and the curva- 
ture of magnetic flux tubes in their small-scale dynamo simu- 
lations; i.e., that the strongest magnetic fields are less curved. 
(This effect could potentially be much stronger in simulations 
similar to ours, but performed at significantly higher spectral 



resolution, where the magnetic diffusion scale would move to 
larger fc-values.) 

We find that the magnetic curvature radius A c and the mag- 
netic rms scale A rm , (cf. Eqs. (16) and (17), respectively, i n|Pa-| 
per 1 1 evolve similarly to A mag . In particular, for the weak-field 
model we find (A c )J'g* « 9.5 km and (A ms )Q'g* « 3.7 km. 
The corresponding values for the strong-field model are about 
10% larger. Note that A c and A rm , combined characterize the 
structure of the magnetic field. They measure respectively 
how sharply magnetic flux tubes are bent and how thinly they 
are stretched. Such information is not contained in A mag alone. 

Spectral kinetic energy distributions from the non-rotating 
models with different initial field strengths are shown in Fig- 
ure QT] The stochastic nature of the SASI and the turbu- 
lent flows necessitates the use of temporally averaged spec- 
tra when cross-comparing the models; in particular, the ki- 
netic energy spectra shown in Figure 1 1 are averaged over the 
time period extending from 800 ms to 1100 ms (that is, we 
plot (e" k in ) o ' 8 s versus fc). The spectra show that the majority 
of the kinetic energy resides on relatively large spatial scales 
(small fc; see the dotted vertical reference lines). The spec- 
tral kinetic energy density is roughly proportional to fc -4 ' 3 
for small fc-values (k < 0.2 km -1 ), while for larger fc-values 
(fc > 0.6 km -1 ) the flow is heavily influenced by numeri- 
cal dissipation and the kinetic energy decreases rapidly with 
increasing wavenumber (e kin ex fc -9 / 2 ). Magnetic field ampli- 
fication is driven by the turbulent flows, and it is the kinetic 
energy of the small-scale motions that is available to be tapped 
by the magnetic fields. 

When comparing the non-rotating models in the left panel 
of Figure 1 1 we see a decreasing trend in the spectral kinetic 
energy density for larger wavenumbers (fc > 0.2) in mod- 
els with a stronger initial magnetic field. (The decrease in 
kinetic energy on small scales is balanced by a correspond- 
ing increase in magnetic energy.) We emphasize this differ- 
ence further in the right panel, where we plot the difference 
in e tin for the stronger field models relative to the weak-field 



16 



ENDEVE ET AL. 



1.0E+00 


■V N 1 . ' 




i iiil 




f iV ^ k- !/J ^" V3 \ 


2Ti/20flx 


2ti/10A* - 












^V \ 








\ ^ ■> 








\ v i 






£ 


\ s ' 






m 






" 










% 1.0E-01 




\ 


-. 


c 


■ V\ 


\ 






>*\ 


\ 


. 










CT- 
















01 








C 




. i 






:rh inni n n nrnn 


'■\s \ 


■ 


<J 




0J 


;B0_lE12LQ_a.OEQO 


%\ 




2 1.0E-02 


= B0 1E13L0 O.OEOO 




v k" V3 - 









Y 








% N \ '■ 












% \ v 


a. 




; A x 






\ V 




i2n/300 km 




1.0E-03 


1 


*■ ■ i 



1.0E-02 



1.0E-01 
k [km"'] 



1.0E + C 



14 


1 !B0_1E12LO. 


I I I | 

.O.OEOO... 




1 1 1 1 1 




B0_1E13LO. 


.O.OEOO 






0.12 


7 






"4 


0.10 


" \ 






- 


0.08 


" \ | 






- 


0.06 


" \ 




I 




0.04 


\ \ 


/ '\ 


/ i l * 


■ ■ 
li 


0.02 


~ \ 


I 
f 

": 1 






onn 


|2tt/J00 km % 


S ^', . 1 


|2Ti/20ax 


2ti/10Ax . 
.... 1 



1.0E-02 



1.0E-01 
k [km"'] 



1.0E + C 



Figure 11. Left panel: (time-averaged) kinetic energy spectra e tm from non-rotating models B10L0.0E00 (solid line), B12L0.0E00 (dotted line), and 
B13L0.0E00 (dotted line). Right panel: difference in spectral kinetic energy relative to the weak-field reference model (1 — e^n/e^), where e k " n f is the 
spectral kinetic energy density of the weak-field model (B10L0.0E00). In both panels we include vertical reference lines to indicate spatial scales of 300 km, 
20 X Al, and 10 X AL In the left panel we also include reference lines proportional to power-laws in k: k~ 5 ' :i (long-dashed), fc~ 4 ' 3 (dash-dotted), and k~ 9 ' 2 
(dash-dot). 



model, 1 — e m /e^. (The spectral kinetic energy density of the 
weak-field reference model is here denoted e^.) The spec- 
tral kinetic energy density in model B12L0.0E00 is reduced 
by up to ~ 6%, while in the strong-field model it is reduced 
by a maximum of ~ 15% relative to the weak-field model. 
The largest difference is seen in the strongly diffusive regime 
around k = 0.6 km -1 . These results demonstrate that even 
relatively weak initial magnetic fields can be amplified and 
impact the flow, although only on small spatial scales. For 
larger spatial scales (k < 0.06) the differences are caused by 
differences in the pre-shock kinetic energy (due to differences 
in R sh and box size L). 

The kinetic energy of the pre-shock flows, density strati- 
fication (resulting in the real-space power-law (in radi us) in 
the spherically averaged kinetic energy density (Section 3.4 1), 
and compressibility (due to the presence of supersonic flows 
below the shock) contribute to the kinetic energy spectra in 
Figure ITT] To further investigate details about the shape of 
the kinetic energy spectra due to these factors, in particular 
the —4/3 slope (as opposed to the —5/3 slope of Kolmogorov 
turbulence), we have computed (1) kinetic energy spectra with 
the pre-shock flow velocity set to zero e^fc), (2) kinetic en- 
ergy spectra with the pre-shock flow velocity set to zero and 
corrected for the radial density stratification e" n (fc) (i.e., we 

use X = \/(p)v L u j with j G {x,y,z} in Eq. (pi), and (3) 
kinetic energy spectra with the pre-shock flow velocity set to 
zero and corrected for radial density stratification, but with 
local compressibility retained ej™(fe) (i.e., we use X = y/puj, 
where p = (p)v L x (p/(p)sVi)i m Eq. (|6|). Results from these 
calculations are plotted in Figure [12J where we plot compen- 
sated kinetic energy spectra from tEe weak-field model (aver- 
aged over the time interval from 804 ms to 918 ms): e kin x fc 5 / 3 
(solid line), e kin x fc 5 / 3 (dotted line), e,™ x fc 5 / 3 (dashed line), 
x /c 5 / 3 (dash-dot line). (In Figure 
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and e™ have 



and e 

been multiplied by a factor of three for convenient comparison 
with e tin and e k ' m .) 
When comparing e 



0.0030 



£ 0.0025 



r, 0.0020 



0.0015 



0.0010 



a- 0.0005 



0.0000 




.0E-02 



.OE-01 
k [km"' 



1.0E+00 



and e kin in Figure 12 we see that the 



Figure 12. Compensated kinetic energy spectra from model B10L0.0E00. 
(The spectra are computed by averaging over the time period from 804 ms to 
918 ms.) The solid line correspon ds to the kinetic energy spectrum shown as 
the solid line in left panel in Figure fTT] The kinetic energy spectrum obtained 
when setting the pre-shock flow velocity to zero e^, is represented by the 
dotted line. Kinetic energy spectra with corrections for density stratification 
(also with the pre-shock flow excluded), e ki ' n and e^", are represented by the 
dashed and dash-dot lines, respectively (see text for details). Note the narrow 
inertial range in the stratification-corrected spectra (Sk';', , e k ™ oc k~ 5 ' 3 ) in 
k £ [0.04, 0.1] km -1 (i.e., spatial scales from ~ 160 km to ~ 60 km). 



supersonic pre-shock accretion flow contributes to the energy 
spectrum, mostly for small fc-values, but the two spectra re- 
main similar in shape. (By integrating the two spectra over all 
k we find E tm w 0.053 B and E^ w 0.040 B.) On the other 
hand, the kinetic energy spectra change markedly when the 
density stratification is excluded from its computation. The 
—4/3 spectral slope seen in the left panel of Figure 1 1 is due 
to density stratification from R ms to R sh . Effects due to com- 
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pressibility are subdominant and do no t change the shape of 
the spectrum in any significant way. ( Kritsuk et al. (2007) 
scaled the velocity with p 1 ' 3 to recover Kolmogorov —5/3 
scaling in spectra from simulations of supersonic isothermal 
turbulence.) Moreover, we observe a narrow inertial range in 
k <G [0.04,0.1] km -1 (i.e., spatial scales from ~ 160 km to 
~ 60 km) where e^, e™ ex fc -5 / 3 . For larger fc-values (around 
k = 0.2 km -1 ) we observe a bump in the e^ and e™ spec- 
tra (e.g., Ishih araet al. |20 09 and references therein), which 
is less pronoun ced for e"„'.~ We note that the simulations by 



Haugen et aI~| ( 2003[J2004| argue in favor of a k~ 
trum for non-helical MHD turbulence. (Given infini 



5/3 



spec- 
(Given infinite spectral 
resolution, the e kin spectrum could possibly also follow —5/3 
scaling for larger k, where the spectrum presumably would be 
less influenced by density stratification.) 

The identification of Kolmogorov-like spectra in our sim- 
ulations helps us associate post-shock flows with turbulence. 
In Figure 12 the peak in the ejj*, and e™ spectra for smaller de- 
values (k w 0.02 km -1 ; i.e., spatial scales around 300 km) 
is associated with the large scale SASI flows (cf. Figure 
[2]), which drive post-shock turbulence. The peak is located 
around k = 0.05 km -1 (i.e., spatial scales around 125 km) 
in the earlier stages of S ASI-development, when the spectrum 
also has a steeper-than — 5/3 slope for larger k. When the 
SASI develops nonlinearly, and the average shock radius be- 
gins to increase, the peak moves to smaller k -values, and the 
e^, and ej™ spectra develop Kolmogorov slopes. Thus, the 
power in the large scale flows cascades to smaller-scale flows. 
In particular, integrating the e^ spectrum in Figure [T2|over k, 
from k = 0.04 to /c max , gives 0.022 B. Thus, a large fraction 
(up to ~ 50%) of the kinetic energy below the shock may be 
associated with turbulence. These observations suggest that 
the SASI saturates due to the development of turbulence via 
secondary instabilities (e.g., the Kelvin-Helmholtz instabil- 
ity), which feed on the power in the low-order SASI modes 



(e.g., Guilet et al. |2010| ). The turbulent energy is either dissi- 
pated via viscous heating, or converted into magnetic energy 
and dissipated via Joule heating. However, we find that sig- 
nificantly less than 50% of the post-shock kinetic energy is 
accessed for magnetic field amplification (cf. Figur e [T3) . 

We use the unmodified kinetic energy spectrum e kin in our 
further analysis since it is related to the total kinetic energy via 
Eq. d9j, and therefore most useful for extracting quantitative 
information from our simulations. 

From the spectral kinetic energy density we obtain the tur- 
bulent kinetic energy in our simulations 



EZ 



(3kin Cbrxi . 



(12) 



For the purpose of studying magnetic field amplification, we 
have chosen to define turbulent flows to include flows resid- 
ing on scales with k > fc, ur = 2n/\ m , where the turbulent 
spatial scale covers 25 grid cells A lul = 25 x Al « 30 km 
(for A/ = 1.17 km). For reference, A lur is more than an order 
of magnitude smaller than the average shock radius, which 
again is comparable to the forcing scale of the turbulent flows 
(i.e., the scale of the supersonic downdrafts from the shock 
triple-point), but comparable to _R PNS (~ 25% smaller). This 
particular choice for fc, ur is motivated by several factors, in- 
cluding (1) most of the magnetic field amplification occurs 
on spatial scales with k > fc, ur (Figures 10 and 14 1, (2) any 
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Figure 13. Turbulent kinetic energy (Eq. \\2) , black lines) and total mag- 
netic energy (red lines) versus time in non-rotating models in which the ini- 
tial magnetic field strength is varied: B10L0.0E00 (solid), B12L0.0E00 
(dashed), and B13L0.0E00 (dotted). The dash-dotted line is proportional 
to exp (£/t), with r = 100 ms. The long-dashed horizontal line indicates 
the upper limit for turbulent kinetic energy, assuming Kolmogorov scaling 
applies for k > k lm . (10~ 2 B = 10 49 erg.) 



k > k m (Figure 111, and (3) the flow Taylor microscale, 



dynamic effect of the magnetic field is seen on scales with 



A T == \/(u 2 )y^/(\V x u| 2 )- ^ h , which m easures the average 
size of turbulent eddies (e.g., Ryu et al. |2000[ ), is comparable 
to A lur (about a factor of two smaller). 

In Figure 13 we plot the time evolution of the turbulent ki- 
netic energy in the non-rotating models. The turbulent kinetic 
energy evolves similarly to the total kinetic energy below the 
shock (cf. top left panel in Figure 15). It grows exponen- 
tially during the ramp-up of the SASI and reaches a satura- 
tion level, where the intermittent time variability is superim- 
posed on a barely noticeable overall growth. During the ex- 
ponential growth phase, the growth rate is somewhat lower 
and the saturation level about an order of magnitude below 
the total kinetic energy beneath the shock. The exponential 
growth time in the weak-field model is r « 100 ms. The 
time-averaged saturation levels for the turbulent kinetic en- 
ergy in the respective models are found to be (■££") o!8* = 
5.90 x 10" 3 B (B10L0.0E00), 5.48 x 10~ 3 B (B12L0.0E00), 
and 5.11 x 10~ 3 B (B13L0.0E00), which implies - 7% 
(4.2 x 10" 4 B) and ~ 13% (7.9 x 10~ 4 B) reductions with 
respect to the weak-field model for model B12L0.0E00 and 
model B13L0.0E00, respectively. These reductions in turbu- 
lent energy are comparable to the increase in magnetic energy 
in the respective models. Thus the magnetic energy grows at 
the expense of the turbulent kinetic energy. 

The saturation level for the turbulent kinetic energy in our 
models is only about a factor of two below what is obtained by 
(hypothetically) assuming Kolmogorov scaling (e ti „ ex fc~ 5 / 3 ) 
for k > fc, ur (indicated by the long-dashed line in Figure [T3). 
Thus, E^ s=s 10~ 2 B may serve as an upper limit for tur- 
bulent kinetic energy in our models, and therefore also as 
a reasonable upper limit on the magnetic energy attainable 
in these simulations. (The turbulent kinetic energy may de- 
pend on the accretion rate ahead of the shock, which is held 
fixed in our simulations. Thus, the upper limit on ££ is only 
approximate.) We find that the magnetic energy in model 
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B13L0.0E00 saturates at about 10% of the turbulent kinetic 
energy. One must keep in mind that the magnetic energy is 
also heavily influenced by numerical dissipation during the 
saturated phase. It is entirely possible that the magnetic en- 
ergy can grow beyond the levels seen in our simulations, but 
probably not much above the long-dashed horizontal line in 
Figure 13 Small-scale dynamo simulations commonly show 



that the magnetic energy spectrum lies slightly above the ki- 
netic energy spectrum on the smallest sca les during saturation 
(e.g., Brandenburg & Subramanian 2005). The magnetic en- 
ergy does not exceed the kinetic energy in any part of the spec- 
trum in our simulations. This may be due to finite resolution 
and numerical dissipation on the smallest scales. 

3.6.2. Varying the Spatial Resolution 



In |Paper I we found that magnetic field amplification from 
SASI-induced turbulent flows is very sensitive to the spatial 
resolution adopted in the numerical simulations. (In gen- 
eral, increased spatial resolution results in stronger magnetic 
fields and improves the conditions for a dynamical influence 
of magnetic fields.) With the energy spectra we continue to 
study the effect of resolution in this section. 

Energy spectra from simulations of the non-rotating model 
with Bo = 1 x 10 12 G for various spatial resolutions are plot- 
ted in Figure [14] In the left panel we plot the spectral kinetic 
energy densityTln the right panel we plot the fractional mag- 
netic energy enclosed by the fc-space shell with radius k 



/«,(*) = J- / e maa (fc')rffc', 

£„, Jk mm 



(13) 



where f„ m (k) is normalized to the total magnetic energy £"„„„, 
so that / mag (fc max ) = 1. Results from three simulations are pre- 
sented, and the grid size has been varied by a factor of two: 
Al = 2.34 km (low resolution; dotted lines), A^ = 1.56 km 
(medium resolution; dashed lines), and A^ = 1.17 km (high 
resolution; solid lines). Kinetic energy spectra are averaged 
over a time period extending from 800 ms to 1100 ms, while 
magnetic energy spectra are plotted for t = 800 ms, 900 ms, 
1000 ms, and 1100 ms. 

The kinetic energy spectra are very similar and follow each 
other closely for small wavenumbers (k < 0.2 km -1 ). Nu- 
merical dissipation influences the kinetic energy for larger de- 
values, and e ki „ begins to fall off more rapidly with increasing 
k. The fall-off starts at smaller fc-values for the lower resolu- 
tion runs: e kin falls below 10~ 3 B km -1 around k = 0.5 km -1 
in the low resolution run, and around k — 0.8 km -1 in the 
high resolution run. 

Since most of the kinetic energy below the shock resides on 
large scales, and includes the flows associated with the su- 
personic stream ahead of the shock triple-point (Figure El, 
the total kinetic energy below the shock, E tiri , is insensitive 
to the spatial resolution. During the highly nonlinear opera- 
tion of the SASI spiral mode we find (_E kin )gg; = 0.045 B, 
0.043 B, and 0.044 B, for the low, medium, and high res- 
olution model, respectively. Our definition of k m , used in 
Eq. (12 1, is not optimal when comparing simulations com- 
puted with different spatial resolutions, since it results in 
smaller /c„ lr and more turbulent kinetic energy in models with 
larger Al ((E™)o.ll = 8.6 x 10~ 3 B, 6.3 x 10~ 3 B, and 
5.5 x 10~ 3 B, for low, medium, and high resolution, respec- 
tively). For the purpose of quantifying the increase in ki- 
netic energy on small scales due to higher resolution, we fix 



fc„ lr = 0.2 km x and find that (-Ev' n )J'8* increases (linearly) 
by a factor of two when A/ is decreased by a corresponding 
factor of two ({EZ}o.ll = 2 -7 x 10~ 3 B, 4.5 x 10~ 3 B, and 
5.5 x 10 -3 B for low, medium, and high resolution, respec- 
tively). 

The self-similar evolution of the magnetic energy spectra is 
clearly displayed in the right panel of Figure [14] The differ- 
ent models are well separated, while each model's temporally 
separated spectra fall practically on top of each other. The 
spectra are shifted to larger /c-values when the resolution is in- 
creased. The shift to the right in the spectrum is a direct con- 
sequence of a corresponding shift of the spatial scale where 
numerical diffusion dominates. The characteristic scale of the 
magnetic field is roughly constant with time over the time pe- 
riod displayed, but decreases linearly with increasing spatial 
resolution. In particular, we find (A mag )g'g' ~ 32 km, 23 km, 
and 18 km, for the low, medium, and high resolution runs, re- 
spectively. Similarly, the magnetic rms scale A mis (the average 
flux tube thickness) decreases by almost a factor of two when 
Al is reduced by a factor of two (from 7.0 km to 3.8 km). 

The shift to smaller spatial scales — in particular the de- 
crease in the flux tube thickness — afforded by higher spatial 
resolution results in stronger magnetic fields and an increase 
in the magnetic energy. The integrated magnetic energy be- 
low the shock in the low resolution model reaches saturation 
for t < 800 ms, and does not grow much beyond 3 x 10 -6 B. 
(Saturation of magnetic energy in this model is solely due to 
numerical dissipation. The influence of magnetic fields on 
small scale flows emphasized in the right panel of Figure 1 1 
is not observed in lower-resolution models.) In the high res- 
olution model the magnetic energy grows throughout the run, 
and reaches 2 x 10~ 4 B near the end (dashed red line in Fig- 
ure [13}, It is interesting to note that during the time span 
from 800 ms to 1 100 ms, between 68% and 77% of the total 
magnetic energy resides on scales smaller than A lur , nearly in- 
dependent of spatial resolution. (There is a weak decrease in 
the percentage with increasing resolution). The correspond- 
ing percentage for spatial scales smaller than 10 x Al varies 
between 16% and 21%. 

We expect the magnetic energy spectra will continue to 
move to higher wavenumbers when the resolution is increased 
beyond that of our simulations. The shift to smaller spatial 
scales (smaller flux tube cross-section) is accompanied by 
stronger magnetic fields, and we expect the flux tube cross- 
section to decrease until the magnetic fields become strong 
enough to caus e a back-reaction on the fluid through the 
Lorentz force. (Haugen et al. (2003) presented converged 
magnetic energy spectra in their simulations of non-helical 
MHD turbulence. In their converged spectra, most of the mag- 
netic energy resides at a wavenumber ~ 5 times the minimum 
wavenumber in the computational domain.) 

3.7. Magnetic Energy Growth Rates 

In this section we focus on the relative importance of mech- 
anisms that control the exponential growth rate of magnetic 
energy when the magnetic field is weak and the kinematic ap- 
proximation remains valid. We also consider the impact of 
finite numerical resolution on the growth rate in our simula- 
tions. 

An eddy turnover time r cildy = A ma „/u|" s is commonly in- 
voked as the characteristic exponential growth time of mag- 
netic fields in a tu rbulent small-scale dynamo (e.g., Kulsrud & 
Anderson 1 1992 1. Here A ma „ is the characteristic spatial scale 
of the magnetic field defined below Eq. ( pT| ). The turbulent 
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Figure 14. Energy spectra from simulations in which the spatial resolution has been varied. Kinetic energy spectra are plotted in the left panel, and the fractional 
magnetic energy enclosed by fc-space shell with radius k, /mag(fc) (cf. Eq. j!3| ), is plotted in the right panel. Results are plotted for the non-rotating model with 
So = 1 X 10 12 (B12L0.0E00). The spatial resolution in these runs has been set to Al = 2.34 km (dotted lines), 1.56 km (dashed lines), and 1.17 km (solid 
lines). The kinetic energy spectra are averaged over the time period extending from 800 ms to 1100 ms. For each model, magnetic energy sp ectra are plotted 
at t = 800 ms, 900 ms, 1000 ms, and 1100 ms (thicker lines represent more advanced time states). The mean magnetic wavenumber (cf. Eq. JTTJ) is denoted 
with a diamond on each spectrum in the right panel. In both panels we include vertical reference lines indicating the spatial scale of 300 km (dash-dot line), and 
10 X Al (with line styles matching each of the models). 



rms velocity is 



2£T 



1/2 



(14) 



where M sh is the mass in V^. The use of M sb in Eq. ( 14 1, in- 
stead of only the mass of the flow included in E^, results 
in an underestimate of u' 1 " . On the other hand, El"' (and 

mis ' kin V 

therefore u'"^) is also sensitive to the definition of k m , which 
may be larger than the value we use in Eq. ( 12 1 and result 
in smaller ££*. Nevertheless, Eq. (14i provides a reason- 
able order-of -magnitude estimate of u™. We find that u™ s 
grows rapidly during the initial ramp up of the SASI and then 
levels off at later times. In the non-rotating models we find 
( u m's)o 8s ~ 4000 km s _1 . A turbulent rms velocity of several 
x 10 3 km s _1 is consistent with an inspection of the subsonic 
flows below the shock: the average velocity among the zones 
with |u|/cg < 1 is about 7000 km s _1 . For the non-rotating 
models r cddy is about 5 ms during the highly nonlinear stage 



of strong SASI activity. (Another commonly used expression 

|2\-l/2 

>V sb ' 



gives a similar 



for the eddy turnover time, (|V x u 
result.) 

We now investigate the individual magnetic energy growth 
rates relevant to our simulations. Assuming a non-ideal elec- 
tric field — (u x B) + r/J with scalar resistivity 77, the evo- 
lution equation for the magnetic energy density is easily ob- 
tained by dotting the Maxwell-Faraday (induction) equation 
with B//xo: 

de 1 

^+V-P = -u-(JxB) B-Vx(TyJ), (15) 

at fj, 

where P = [u(B-B)-B(B-u)]//i and J = (V x B)/fj, . 
(See als o Eq. (10) as well as the discussion in Section 3.3 in 
Paper I ) The first and second terms on the right-hand-side of 
Eq. (T3 1 represent work done against the Lorentz force (W L ) 
and magnetic energy decay due to resistive (Joule) dissipation 



(— Qj), respectively. Kinetic energy of the flow is converted 
into magnetic energy if W L > 0. 
It is apparent from Eq. (fl~5]> that the total magnetic en- 



ergy growth rate t„ = 



v {de imi /dt) equals the sum 



'JxB~ 



-t s l of individual rates due to work done against 
the Lorentz force, accretion of magnetic energy (Poynting 
flux) through dV ms , and resistive energy dissipation. (The an- 
gle brackets in the total rate imply an integral over a volume 
bounded by the surface dV Pm and a spherical surface enclos- 
ing the accretion shock.) The Poynting flux through the spher- 
ical surface enclosing the accretion shock vanishes because 
u || B ahead of the shock . The Poynting flux through dV PNS 
and resistive dissipation generally result in decay of the mag- 
netic energy in the computational domain. The decay must be 
overcome by the Lorentz work term in order for the magnetic 
energy to increase. 



In Paper I we found flux tube stretching by turbulent flows 
driven by the spiral SASI mode to be the dominant mechanism 
for magnetic field amplification (see also Figure 15 below). 
The magnetic energy growth rate due to work done against 
the Lorentz force is 



'JxB 



— / u • (J x B) dV 



-E-mug JV 

: 2< s /A milg = 2r 



-1 

*idy 1 



(16) 



where the turbulent rms velocity u l ™ s and the characteristic 
scale of the magnetic field A m .,„ have been used. (The factor 



two in the second part of Eq. ( 16 1 stems from the factor of one 
half in the definition of magnetic energy, but is probably not 
important for this rough estimate.) The corresponding growth 
time is then approximately 



tjxb ~ 2.5 ms 



A™ 



20 km 



4000 km s~ 



(17) 
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Figure 15. Magnetic energy growth rates versus time for the non-rotating 
weak-field mo del ( BIOLO.OEOO). The growth rates are based on terms ap- 
pearing in Eq. J15) and are due to work done against the Lorentz force (solid 
line) and Poynting flux losses due to accretion through the spherical surface 
with r = Rpns (dash-dot line). We also plot the growth rates due to compres- 
sion (dotted line) and stretching (dashed line) (Eqs. (13) and (1 1) in Paper I 
respectively). 



The volume occupied by the PNS is excluded from our sim- 
ulations, and the magnetic energy in V sh is also affected by 
accretion of magnetized matter through dV ms . The decay rate 
due to this process is 



1 



P-dS 



0\h 



3M 

2tt p„L 3 B 



(18) 



where in the rightmost estimate we have adopted the expo- 
nential decrease of magnetic field with radius over a charac- 
teristic length scale Lb (cf. Figure [6J to relate the magnetic 
energy in V^ h to the field strength at the surface of the PNS: 

E m 

the inner boundary is then approximately 



Y^^-Lg. The decay time due to accretion through 



rp « 90 ms 



Po 



3 x 10 10 g cm 



-3 



( Li 



M 



V 100 km/ \O.36M s- 



(19) 



The average mass density around r = R ms , denoted po, stays 
fairly constant throughout the simulations. 
In Figure 15 we plot the growth rates rJ^ B (solid line) and 



Tp 1 (dash-dot line) versus time for model B10L0.0E00. This 
model exhibits exponential magnetic energy growth through- 
out with a growth time of about 66 ms. The growth rates 
plotted in Figure [15] are computed from numerical approxi- 
mations (second-order finite differences) to the integral ex- 
pressions, and not the approximations provided by the right- 
most expressions in Eqs. (17) and (19). We also include 



the growth rates due to stretching Ty u and compression Ty. u 
(Eqs. (11) and (13) in Paper I, respectively), and the plot 
shows that stretching dominates over compression. The rates 
remain quasi-steady for t > 750 ms, and in particular we 



find 



ii.i. 

'JxB/0.9. 



480 



and (t, 



-lu.i, 

• /0.9s 



9 



(We 



also find 



'Vu 



l.i. 

'0.9* 



515 s" 



and 



a.i, 

/0.9s 



76 



•) 



We note that there is good agreement between the numeri- 
cally computed growth rates and the growth rates predicted 
by the estimates provided by the rightmost expressions in Eqs. 
( |17[ ) and ( 19 ). Furthermore, the relative importance of these 
rates in determining the total magnetic energy growth rate be- 
comes clear: since tj x b <C Tp, accretion of magnetic energy 
through dV PNS has virtually no effect on the growth of mag- 
netic energy in V sh . 

The discrep ancy between the millisecond growth time pre- 
dicted by Eq. ( 16 1 and the numerically measured growth time 
(t sa 66 ms; Figure [5) suggests that numerical dissipation 
plays an important role in controlling the growth time for 
magnetic energy in our simulations. Thi s is further supported 
by the results presented in Section 3.6 which show that the 
magnetic energy develops on spatial scales that are strongly 
affected by numerical dissipation (see also Appendix [A}. If 
not suppressing field growth entirely, numerical dissipation 
tends to increase the magnetic energy growth time. The char- 
acteristic decay rate due to resistive dissipation of magnetic 
energy is 



1 

El 



1 
v A*o 



B • V x (77J) dV 



2q 
A?' 



(20) 



where we have introduced the dissipation scale A d . The decay 
time due to resistive dissipation is then 



R„, 



A mM 



r JxB, 



(21) 



where the magnetic Reynolds number is defined as R m = 
"l^/l- The magnetic Reynolds number in the supernova 
environment is expected to be extremely large ; on the order 
10 17 in the PNS ( |Thompson & Duncan Il993) . As far as the 
magnetic energy growth rate is concerned, resistive effects are 
only relevant on very small scales, and the growth is most 
likely curbed by dynamical interactions with the fluid through 
magnetic tension forces before the magnetic field develops 
to resistive scales du e to turb ulent stretching of flux ropes 
(Thompson & Duncan 1993). Thompson & Duncan (1993) 
list (in their Table 1) the resistivity in the PNS convection zone 
(77 = 1 X 10~ 4 cm 2 s _1 ). Adopting this value, the resistive 
decay time for a magnetic field varying on a spatial scale of, 
for example, 1 m (i.e., much smaller than any scale resolved 
by our simulations) becomes very long (7) = 5 X 10 7 s) com- 
pared to the explosion time (~ 1 s). 

Resistive effects are, however, important to consider in nu- 
merical MHD simulations of astrophysical systems. We do 
not explicitly include resistivity in our simulations, but the 
numerical scheme incorporates an effective numerical resis- 
tivity in the induction equation in order to stabilize the solu- 
tion when discontinuities or underresolved gradients appear 
in the flow (see Appendix [A] for further details). An approxi- 
mation of the total growth rate can be obtained by combining 
Eqs. ([16]) and (|20|: 



'JxB 



1 - 



1 (\ 



Rr, 



A, 



(22) 



In our simulations we have r lu , 3> tj x b and A d < A mag . De- 
fined this way, the magnetic Reynolds number in our simula- 
tions is therefore somewhat larger than, but still close to, unity 
(R m > 1). This conclusion is consistent with the observations 
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from the energy spectrum plots above, which show that a siz- 
able fraction of the magnetic energy resides on spatial scales 
where numerical diffusion is significant. Our simulations are 
therefore likely to grossly underestimate the magnetic energy 
growth rates that can be expected under more realistic physi- 
cal conditions (i.e., where R m 3> 1). 

We point out here that the saturation of magnetic energy 
observed in model B13L0.0E00 does not mean that rJ" xB ~ 
s _1 for this model. We find that the amplified magnetic 
fields in model B13L0.0E00 result in about a 10% reduction 
in Tj^g relative to the weak-field model, which, because of 
numerical dissipation, results in a significant reduction in the 
total growth rate r" 1 (we expect t, -1 to be the same in both 
models). 

In Appendix [A] we measure numerically the magnetic en- 
ergy decay rate due to resistive dissipation in one of our sim- 
ulations. We find r~ l « 380 s _1 , which is comparable to, 
but still somewhat smaller than tJ^b ■ The decay rates t~ 1 , 
along with A d (« A mag ) and rj mm (Appendix Ab, and tJ^-q do 
not vary significantly with time during the nighly nonlinear 
stage of the SASI. Thus, the numerically measured growth 
time (r" 1 s=s 66 ms; Figurepjl is mostly the result of two large 
and competing processes: growth due to tJ^-q and decay due 
to T~ l , For increasing spatial resolution (i.e., increasing R m ) 
we expect r" 1 — > Tj^ B (cf. Eq. (|22Jl). 

Indeed, we have carried out simulations with different spa- 
tial resolutions and measured the magnetic energy growth rate 
when an epoch of exponential growth can be identified. We 
find that the growth rate increases with increasing resolution: 
for the lowest resolution model (Al — 2.34 km) the exponen- 
tial growth time is about 150 ms, while in a simulation with 
Al = 0.78 km the exponential growth time decreases to about 
50 mrj In fact, a divergent increase in the magnetic energy 
growthrate with increasing magnetic Reynolds number (i.e., 
resolution) has been reported in direct n umerical simulations 
of MHD turbulence (e.g., Haugen et al. 2004) and recently in 
simulations of turbulent star formation using adaptive mesh 
refinement (e.g., Syr et al. 2010 Federrath et al. |2011|l. 



These authors show results from simulations in which the res 
olution has been doubled several times, and they find that the 
magnetic energy growth rate increases as a power law with 
increasing magnetic Reynolds number. 

The sensitivity of the magnetic field evolution to numerical 
resolution does raise concerns about what aspects of our sim- 
ulations are relevant to core-collapse supernovae. Dissipation 
due to finite grid resolution tends to suppress magnetic energy 
growth. At face value, our simulations (falsely) result a neg- 
ative assessment on the efficiency of SASI-induced magnetic 
field amplification. However, as the resolution is increased 
the growth rate increases and the resulting magnetic fields be- 
come stronger. Further analysis suggests that the simulations 
grossly underestimate the growth rates and fields that may ob- 
tain in the supernova environment. The SASI-induced turbu- 
lent magnetic field amplification mechanism is a robust result 
from our simulations. Only the growth rate, saturation am- 
plitude, and the dynamical impact of the amplified magnetic 
fields remain uncertain. An important consequence of the im- 



To conserve computational resources, this model was not run to com- 
pletion but until the computational domain consisted of 1280 3 zones (t Ri 
880 ms). At this time the SASI is still ramping up and the magnetic energy 
growing rapidly. 



plied millisecond growth time is that any weak seed magnetic 
fields may be amplified to saturation levels (|B| w y//iQp\u\) 
in a core-collapse supernova if the SASI operates and drives 
vigorous turbulent flows below the shock. The kinetic energy 
available to amplify the magnetic energy (some fraction of 
££*) is not sufficient for magnetic fields generated in this way 
to play a principal role in the explosion dynamics. We can- 
not, however, completely rule out the possibility that SASI- 
generated magnetic fields play a secondary role in the dynam- 
ics leading to core-collapse supernovae. 

3.8. Simulations with Initial Rotation 



Results from rotating models are shown in Figure 16 in 
which we plot kinetic energy (left panel) and magnetic en- 
ergy (right panel) versus time. Rotating models with Bo = 
1 x 10 10 G, and l = 1.5 x 10 15 cm 2 s" 1 (B10L1.5E15; 
dashed lines) and l = 4.0 x 10 15 cm 2 s" 1 (B10L4.0E15; dot- 
ted lines) are compared with the corresponding non-rotating 
model (B10L0.0E00; solid lines). 

The most notable difference between these models is the 
earlier onset of the SASI observed in the rotating models. The 
post-shock flow is set into rotation about the z-axis as the 
pre-shock material with angular momentum advects down- 
stream. The kinetic energy in the post-shock flow in model 
B10L1.5E15 increases initially by ~ 50%, and settles mo- 
mentarily into a quiescent state, which lasts for about 200 ms. 
Then, for t > 300 ms, the nonlinear phase of the SASI sets 
in, and the kinetic energy begins to grow exponentially with 
a growth time r w 55 ms, which is notably faster than in the 
non-rotating model. Model B10L4.0E15 receives a stronger 
initial perturbation due to more angular momentum ahead of 
the shock, and the kinetic energy in this model grows rapidly 
by a factor of ~ 5 before settling into a short, quasi-steady 
state with E ki „ ~ 10 49 erg. The kinetic energy begins to grow 
again for t > 200 ms. The kinetic energy in all models even- 
tually reaches similar levels in the strongly nonlinear phase 
of the SASI; when averaged over the last 200 ms of each run 
we find (£L)8-§1; = 0.052 B and (.EL) oil! = °- 050 Bs for 
models B10L1.5E15 and B10L4.0E15 respectively. (We re- 
0.051 B for model B10L0.0E00 in Section 



,1.1. 



jorted (E, 

Hi) 

The earlier onset of the nonli near SASI in th e rotating mod 
els is consistent with Blondin & Mezzacappa ( 2007 1 and Ya- 



masaki & Foglizzo (2008). However, model B10L4.0E15ls 
perturbed relatively hard when the rotating pre-shock material 
advects downstream, and the model does not settle into a qui- 
escent state, as is observed in B10L1.5E15. We think it is very 
likely that the early SASI development in model B10L4.0E15 
is partially a result of our method of initiating the rotating 
models. Nevertheless, the purpose of these simulations is to 
study the effect of rotation on turbulent magnetic field ampli- 
fication during the non-linear phase, and our rotating models 
are suitable for this purpose. 

The evolution of the magnetic energy below the shock dur- 
ing nonlinear SASI-operation (right panel in Figure 16 1 in 



the rotating models is not significantly different from model 
B10L0.0E00. All models exhibit exponential magnetic en- 
ergy growth during the late stages. The magnetic energy in 
model B10L1.5E15 grows exponentially with a growth time 
r w 44 ms during the early stages (from t sa 340 ms to t sa 
550 ms), and grows at a rate similar to model B10L0.0E00 
later on (t > 600 ms). The magnetic energy in model 
B10L4.0E15 grows exponentially at a somewhat slower rate 
than the other models (r « 85 ms). However, all models 
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Figure 16. Kinetic energy (left) and magnetic energy (right) versus time from simulations with varying degree of initial rotation. The specific angular momentum 



in the pre-shock flow has been set to I = 0.0 cm 2 s" 1 (solid; BIOLO.OEOO), 1.5 X 10 15 cm 2 s" 1 (dashed; B10L1.5E15), and 4.0 X 10 15 . 
B10L4.0E15). The magnetic field strength at r = i? PN s is initially Bo = 1 X 10 10 G in all the models. 

have reached similar levels at the end of the respective runs. 
In particular, we find E m , e « 1.2 x 1(T 7 B (t = 878 ms) 
for B10L1.5E15 and E„ m « 3.4 x 1(T 8 B (t = 678 ms) for 
B10L4.0E15. 

4. MAGNETIZATION OF PROTONEUTRON STARS 



(dotted; 



In |Paper I| we pointed out that the underlying PNS may be 
significantly magnetized due to SASTinduced magnetic field 
amplification. In this section we estimate in a similar manner 
the degree of PNS magnetization predicted by the current set 
of simulations. 



Model 


Table 2 

PNS Magnetic field estimates 

t Md (ms) B mag ,p (Erg) 


(B PNS ,P> (G) 


BIOLO.OEOO 
B10L1.5E15 
B10L4.0E15 


1100 
878 
678 


1.14 x 10 44 
6.70 x 10 43 
2.52 x 10 43 


3.3 x 10 12 
2.5 x 10 12 
1.5 x 10 12 


B12L0.0E00 
B12L1.5E15 
B12L4.0E15 


1126 
1000 
644 


3.16 x 10 4/ 
1.15 x 10 48 
1.74 x 10 47 


1.7 x 10 14 
3.3 x 10 14 
1.3 x 10 14 


B13L0.0E00 


1100 


4.49 x 10 4S 


6.5 x 10 14 



Adopting Eq. ( 15 1, the magnetic energy in the volume oc- 
cupied by the PNSTp NS at some time t > to is 

£ m ., g (i) = £, mg (M 

dt' (f (W L - Q,) dV 

1 to \J VpNS 



Note. — Magnetic energy accumulated on the proto-neutron star in com- 
puted models. The inferred magnetic field strength resulting from SASI- 
induced magnetic field amplification is also listed. 



P-dS 



i)\],s 



(23) 



Here P is the Poynting flux through the surface of the PNS, 
and W L and Q s are obtained from u and B, which must be 
computed with an appropriate physical model of the PNS. 
The dissipative term Q s also involves the resistivity ij. Re- 
sistive dissipati on is not likely to su ppress field amplification 
in the PNS (Thompson & Duncan 1993), but may be impor- 
tant to the long-term evolution of the neutron star magnetic 
field (strength and topology). 

Evaluation of the volume integral on the right-hand-side 
of Eq. (23 i involves numerical simulations of the hydro- 



2003); and convective dynamo action, driven by entropy gra- 
dients lepton gradients, or both ( e.g., [Thompson & Duncan 
l"99"3"||Bonanno et al. |2003][2005] >. All these mechanisms op- 
erate inside or on the surface of the PNS and rely on rotation. 
We exclude the PNS from our simulations and do not address 
field amplification mechanisms in its interior. Our simula- 
tions, however, focus on field amplification by the SASI exte- 
rior to the PNS, which is often ignored in models addressing 
the origin of pulsar magnetism. 

From our simulations we compute the increase in magnetic 
energy in the volume occupied by the PNS due to the Poynting 
flux through the surface bounding it, 



£-.,p(i) 



dt' 



P-dS. 



(24) 



m pv 



magnetic evolution inside the PNS during the explosion phase 
of core-collapse supernovae and subsequent PNS cooling, 
and includes neutrino radiation-magnetohydrodynamic sim- 
ulations of dense nuclear matter. Such calculations are well 
beyond the scope of this study. Earlier works have suggested 
numerous mechanisms for field amplification in the PNS, in- 
cluding winding by differential rotation (e.g., Wheeler et al. 
2002); the magneto-rotational instability (Akiyama et al. 



We then estimate the PNS magnetic field due to SA SI activity 
(Bpns,p) = (2^E m .^ P /V ms ) 1 / 2 (cf. Eq. (18) in |Paper l|. 



Results from these estimates for the rotating and non-rotating 
models with varying initial magnetic field strengths are listed 
in Table [2] 

Our results show that the magnetic energy generated by 
SASI activity may result in significant magnetization of the 
PNS. The magnetic energies generated in some of the models 
meet the energy requirements to power the total flare energy 
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released per SG R and the persistent X-ray emission (Thomp- 
|son & Dun can 2001). The models with the weakest initial 
magnetic field predict field strengths in the range of ordinary 
pulsars (a few x 10 12 G), while the models with stronger ini- 
tial magnetic fields predict fields in the magnetar range (ex- 
ceeding 10 14 G). The magnetic field in the strong-field model 
(B13L0.0E00) saturates dynamically, and this model may rep- 
resent an upper limit to the fields attainable from this process. 
On the other hand, the weak-field models do not reach sat- 
uration. The magnetic energy in these models continues to 
grow at an underestimated rate, and the maximum attainable 
field strength/energy is also limited by finite grid resolution. 
The PNS field strengths predicted by these models are there- 
fore artificially low. It then seems likely, given infinite grid 
resolution, that PNS magnetic fields can exceed 10 14 G due 
to the SASI alone, independent of the initial magnetic field 
strength. Moreover, as finite grid resolution severely limits 
the exponential growth rate of magnetic energy, the duration 
of SASI operation may be less critical. The amount of initial 
rotation in the models does not seem to affect the degree of 
PNS magnetization. 

The field strengths listed in Table|2]should also be corrected 
for additional magnetic field amplification as the PNS cools 
and contracts. From conservation of magnetic flux through 
the PNS surface, contraction from a 40 km radius to a radius 
of about 15 km boosts the surface field by a factor of ~ 7. 

We point out that the PNS magnetic fields resulting from the 
turbulent flows driven by the SASI are likely small-scale and 
disordered. A connection to the dipolar magnetic field struc- 
ture inferred for neutron stars is currently missing, and, of 
course, the SASI alone cannot fully explain the origin of pul- 
sar magnetism. However, pulsar magnetic fields are thought 
to consist of a global dipole field superimposed with higher 
order multipole (small-scale) components, and pulsar mag- 
netism is most likely a result of the combined action of mul- 
tiple amplification mechanisms. While the inferences we can 
make are limited by resolution (which affects the magnetic 
growth rate) and the absence of important physics (which de- 
termines the time to explosion), our simulations suggest that 
the SASI could in principle make a nontrivial contribution. 

5. SUMMARY, DISCUSSION, AND CONCLUSIONS 

We present results from three-dimensional MHD simula- 
tions of the SASI. The simulations are initiated from a con- 
figuration that resembles the early stalled shock phase in a 
core-collapse supernova, albeit with simplified physics that 
excludes critical components of a supernova model (e.g., neu- 
trino transport, self-gravity, and the PNS itself). On the other 
hand, our simulations are computed with a spatial resolu- 
tion that is currently inaccessible to state-of-the-art supernova 
models in three spatial dimensions, and they may therefore 
provide valuable insight into MHD developments in core- 
collapse supernovae. In particular we study the evolution and 
amplification of magnetic fields in S ASI-driven flows in order 
to assess the effects of the amplified magnetic fields on super- 
nova dynamics, and the possible role of the SASI in magne- 
tizing the PNS. 

This paper is a continuation and extension of the study initi- 
ated in Paper I The simulations reported here were performed 
with higher spatial resolution (up to 1280 3 grid cells) and 
cover a broader parameter range than the 3D simulations pre- 
sented in Paper I we have varied the strength of the initial 
magnetic field and the degree of rotation in the flow ahead 
of the shock (including no rotation). We have also varied the 



spatial resolution in some of the simulations, and extended 
the analysis from Paper I to include a Fourier decomposition 
of the kinetic energy and magnetic energy in the simulations. 
Our main findings are 

1. The S ASI-driven turbulence that develops is essen- 
tially non-helical, and shares similariti es with convec- 
tively drive n MHD turbulence (e.g., Brandenburg et 



al. 



199 6]). (See also "box turbulence" simulations 
by Haugen et al. ( |2004 i.) When corrected for den- 
sity stratification, the kinetic energy spectra associ- 
ated with the post-shock flow develop Kolmogorov- 
like —5/3 scaling (i.e., ej^, e™ oc fc~ 5 / 3 ; Section 3.6 1 
in a narrow wavenumber range. Moreover, inspection 
of the time evolution of the kinetic energy spectra re- 
veals that the power in low-order SASI modes (i.e., 
large scale flows) cascades to higher-order modes (i.e., 
smaller-scale flows), and that a significant fraction (up 
to ~ 50%) of the post-shock kinetic energy can be as- 
sociated with turbulence (although a smaller fraction is 
involved in magnetic field amplification). This further 
suggests that the non-linear SASI saturates due to the 
development of turbulence via secondary instabilities 
(e.g., |Guiletetal. |2010| i. 

2. The magnetic energy grows exponentially with time in 
turbulent flows driven by the SASI, as long as the kine- 
matic regime obtains. Our simulations de velop flows 
cha racteristic of the SASI spiral mode (e.g., |Blondin &| 
Mezzacappa 2007). These flows drive vigorous turbu- 
lence below the shock (u l ™ s ~ 4000 km s _1 ), which am- 
plifies magnetic fields by stretching. The resulting mag- 
netic field is highly intermittent and consists of thin, in- 
tense magnetic flux ropes. 

3. Simulations initiated with weak or moderate rotation 
evolve similarly to non-rotating models as far as the 
magnetic field amplification mechanism is concerned. 
However, models with initial rotation develop the non- 
linear spiral SASI flows earlier, and exponential mag- 
netic energy growth sets in sooner. The earlier onset 
of the SASI in models with initial rotation is consistent 



with the re sults of Blondin & Mezzacappa (2007) and 
Yamasaki & Foglizzo (2008). 



The magnetic energy grows at the expense of the kinetic 
energy available in the turbulent flows driven by the 
SASI. Our simulations show that strong magnetic fields 
emerge on small (turbulent) spatial scales, and reduce 
the turbulent kinetic energy on those scales. For our ref- 
erence spatial resolution, magnetic fields impact flows 
on scales with wavenumber k > fc dy „ w 0.1 — 0.2 km^ 1 
(A dyn = 27r/fc dyn < 30 — 60 km) and peak around k = 



0.6 km -1 (~ 10 km) (Figure [111). That is, magnetic 
fields do not affect the portion of the kinetic energy 
spectrum with k < /c dy „. The turbulent kinetic energy 
(that is, the kinetic energy on spatial scales below some 
specified cutoff) in models with larger magnetic fields 
is reduced compared to models initiated with weaker 
magnetic fields, indicating a dynamical impact of the 
amplified magnetic field. 

5. The magnetic field evolution in our simulations remains 
very sensitive to the spatial resolution. Key parameters 
extracted from simulations performed with increasing 
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spatial resolution do not converge in the range covered 
in this study. Both the final magnetic energy attained 
and the rate at which the magnetic energy grows in- 
crease with increasing grid resolution. In particular, es- 
timates using data extracted from our simulations sug- 
gest that the magnetic energy may grow exponentially 
on a millisecond timescale under physically realistic 
conditions, with very large magnetic Reynolds num- 
bers, as opposed to the ~ 50-60 ms timescale measured 
directly in our runs. 

The magnetic energy saturates when the magnetic en- 
ergy density becomes comparable to the kinetic en- 



ergy density (i.e., 



|u|) in localized re- 



gions of the flow. Only our "strong-field" model (with 
the largest initial magnetic field) reaches this saturated 
state. The subsequent magnetic field evolution remains 
highly dynamic: strong fields are advected through the 
flow, are temporarily weakened, and then reemerge in a 
seemingly stochastic manner. 

7. The magnetic fields amplified by the SASI are not 
likely to play an important role in the explosion dy- 
namics (but see further discussion below). The pres- 
ence of amplified magnetic fields does not result in no- 
ticeable effects on the global shock dynamics in our 
simulations, and this can be understood as a matter of 
simple energetics. Magnetic energy grows at the ex- 
pense of kinetic energy, and the kinetic energy content 
in the post-shock flow during vigorous SASI activity 
(~5x 10~ 2 B) is not enough for magnetic fields to be- 
come energetically significant to the explosion (~ 1 B). 
This was also pointed out in Paper I Moreover, the 
turbulent kinetic energy — which powers SASTdriven 
field amplification — accessible for magnetic field am- 
plification only amounts to about 10% of the total ki- 
netic energy below the shock. We further point out that 
our estimate for turbulent kinetic energy is not criti- 



cally sensitive to the numerical resolution (Section 3.6 
and Figure [T3). A rapidly rotating (millisecond period) 



PNS would provide an energy reservoir large enough to 
power ma gnetically-driven explosions (e.g., Burrows et 
al. |2007 ), but it is not likely that rotation would be this 



strong in most supernova progenitors (e.g., Heger et al 
2005). These observations suggest a rather passive role 
of magnetic fields in the overall dynamics of at least 
most supernovae. 

8. Our simulations suggest that SASTinduced magnetic 
field amplification may play an important role in de- 
termining the strength of the magnetic field in proto- 
neutron stars and young pulsars. Upon integrating 
the Poynting flux through the surface encompassing 
the PNS, we estimate that the magnetic energy accu- 
mulated on the PNS may account for magnetic field 
strengths exceeding 10 14 G. This is stronger than the 
canonical dipole field inferred for typical pulsars, and in 
this connection two points must be emphasized. First, 
SASTdriven amplification is expected to cease when 
the explosion takes off, so that different delay times to 
explosion (which may for example be a function of pro- 
genitor mass) may result in different degrees of PNS 
magnetization. Second, the SASTamplified portion of 
the field accumulated by the PNS will at least initially 



be disordered and not of the large-scale, dipolar char- 
acter of the fields inferred from pulsar spindown. 

Despite the pessimism of point 7 above regarding the rel- 
evance of SASTamplified magnetic fields to the explosion 
mechanism, we caution that the sensitivity of magnetic field 
amplification and evolution to numerical resolution prevents 
us from completely dismissing magnetic fields as unimportant 
to supernova dynamics in weakly rotating progenitors. 

Certainly, our simulations cannot accurately describe the 
dynamical interaction between the magnetic field and the fluid 
on small scales. An initially weak magnetic field is ampli- 
fied exponentially in turbulent flows when the flux tubes are 
stretched and their cross sectional area decreases. In a realis- 
tic post-shock supernova environment, where R m is extremely 
large, field amplification is likely quenched by dynamic back- 
reaction on the fluid before the flux tube thickness reaches 



the resistive scale (Thompson & Duncan 1993). The resis- 
tive decay time then remains much longer than the dynamical 
timescale of hydro-magnetic interactions. But in numerical 
simulations the flux tube cross section inevitably approaches 
the grid scale, and numerical dissipation sets in and prevents 
further strengthening of the magnetic field. This occurs in 
all our simulations. (The strong-field model (B13L0.0E00) 
develops dynamically relevant magnetic fields, but is also 
strongly affected by numerical dissipation.) 

Our simulations suggest that magnetic fields become dy- 
namically relevant on spatial scales smaller than A dy „ ~ 30 km 
(Figure [lT). The global shock dynamics remains unaffected 
by the presence of magnetic fields (e.g., Figure 15). However, 
we cannot rule out the possibility that flows on scales larger 
than X dyn could ultimately be affected by hydro-magnetic in- 
teractions emerging from small-scale turbule nt flows. Simu- 
lations of non-helical MHD turbulence (e.g., |Haugen et al. | 
2004 ) show that the magnetic energy grows exponentially on 
the turnover time, on all spatial scales during the kinematic 
regime. (We also observe exponential growth on all scales 
in our runs during this regime.) The kinematic regime ends 
when the magnetic energy becomes comparable to the kinetic 
energy. This occurs on a scale by scale basis. Magnetic en- 
ergy growth slows down considerably after this equipartition, 
which occurs first on the smallest spatial scales, and the mag- 
netic energy spectrum settles somewhat above the kinetic en- 
ergy spectrum. (We also observe that magnetic energy growth 
is quenched when e, mg ~ e kin , but the magnetic energy spec- 
trum stays below the kinetic energy spectrum for all k in our 
simulations.) At later times in MHD turbulence simulations, 
the largest spatial scale at which e ma „ > e tin (i.e., A dy „) in- 
creases, and may approach the driving scale of the turbulent 
forcing. For helical MHD turbulence, which may be more rel- 
evant when a rapidly rotating PNS is included in the model, 
A dyn can even grow beyond the forcing scale (e.g., Meneguzzi 
et al. 1981; Brandenburg 2001). However, the timescale for 
this process is relatively slow, and increases with R m (Bran- 
denburg |2001|l. Nevertheless, it would be desirable to deter- 



mine the largest scale at which the magnetic energy equili- 
brates with the kinetic energy in SASTdriven flows. The lack 
of sufficient spectral resolution in our simulations prevents us 
from determining whether magnetic fields can become strong 
on large enough scales to alter post-shock flows in a signifi- 
cant way. 

The SASI may play an important role in improving the 
conditions for successful neutrino -drive n explosion s (e.g., 
|Bruenn et al. ||2006[ |Buras et aT~||2006| pvlezzacappa et al.| 
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112007] |Marek & Janka |[2009l |Suwa et aT~~||20T0l |Mtffleretl 
al. |2012 i. If amplified magnetic fields can alter the evolution 
of the SASI and change the conditions (making them more, 
or less, favorable) for energy deposition by neutrinos, then 
magnetic fields may play a secondary but relevant role in the 
dynamics of a broader range of core-collapse supernovae (i.e. 
not just those arising from rapidly rotating progenitor stars). 



This point was also argued by Obergaulinger & Janka ( 2011) , 
who studied magnetic field amplification in non-rotating col- 
lapsed stellar cores using axisymmetric simulations that in- 
cluded the PNS and neutrino transport. They found that the 
SASI and convection contribute to magnetic field amplifica- 
tion, and observed the most pronounced shock expansion in 
the model where the magnetic field was strong enough to al- 
ter the post-shock flow topology. (This model was initiated 
with a strong pre-collapse magnetic field.) 

However, axial symmetry severely constrain s magnetic 
field evolution driven by the SASI (see Paper I) and, most 
likely, also convectivel y driven fiel d amplification. Simula- 
tions similar to those of Obergaulinger & Janka ( 201 1 1 in full 
3D, where the SASI spiral mode can develop and drive tur- 
bulent field amplification, are therefore highly desired. Such 
simulations will improve on our simulations in (at least) two 
important ways: 

1 . A significant amount of magnetic energy (comparable 
to that in V sh ) is lost through the boundary at r = R PNS 
in our models, and not accounted for in the subsequent 
dynamics. 3D simulations with the PNS included do 
not suffer from this artificial limitation, and will allow 
us to better assess the role of SASI-induced magnetic 
fields. 

2. Simulations that include neutrino transport develop 
neutrino-driven convection, both in the PNS and in the 
shocked mantle. This convective activity will impact 
the evolution of magnetic fields, and possibly also the 
SASI. We will then be able to study magnetic field evo- 
lution in a much more physically realistic supernova en- 
vironment. Moreover, with neutrino transport included, 
we will be able to directly address the role of magnetic 
fields on neutrino-powered explosions. 

The constraint on numerical resolution in order to prop- 
erly describe turbulent flows may still be computationally 
prohibitive, especially when additional (necessary) physics 
components are added to the models. This may be par- 
tially circumvented with the use of adaptive mesh refinement 
techniques and improved numerical algorit hms. Local (o r 
semi-global) simulations (e.g., Obergaulinger et al. 2009), 
adopting physical conditions and forcing functions relevant 
to the supernova environment (i.e., derived from global multi- 
physics simulations), may also be necessary to study turbulent 
magnetic field evolution and its impact on supernova dynam- 
ics in more detail. More investigations, using both local and 
global simulations, are needed to better understand the role of 
magnetic fields in core-collapse supernovae. 

In summary, we conclude from our simulations that mag- 
netic fields in core-collapse supernovae may be amplified ex- 
ponentially by turbulence on a millisecond timescale; i.e., 
much shorter than the time between core bounce/shock for- 
mation and initiation of the explosion. Details of the impact 
on explosion dynamics by SASI-amplified magnetic fields re- 
main unclear, but on energetic grounds alone the role of these 



magnetic fields is likely sub-dominant. The simulations fur- 
ther suggest that small-scale neutron star magnetic fields in 
the 10 14 — 10 15 G range may be formed, which may be suffi- 
cient to power some of the energetic activity that define AXPs 
and SGRs. 
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APPENDIX 

DISSIPATION OF MAGNETIC ENERGY IN THE NUMERICAL SIMULATIONS 

Here we briefly describe the dominant source of magnetic energy dissipation in our numerical simulations. We evolve the 
MHD equations with a second-order, semi-discrete, central-u pwind, finite volume scheme for hyperbolic conservation laws, 
combined with the constrained tran sport (CT) method of Evans & Hawley ( 1988) for divergence -free magnetic field evolution 
(see Kurganov et al. 2001, Londrillo & Del Zanna 2004 and the references there in; in particular, the MC-HLL-UCT scheme in 
Londrillo & Del Zanna ( 2004)). Furthermore, we adopt the HLL Riemann solver ( Hart enet al. |1983 i to compute the fluxes and 
electric fields needed to evolve the system of MHD equations. The HLL Riemann solver considers only the fastest left- and right- 
propagating characteristic waves of the underlying hyperbolic system (fast magnetosonic waves for MHD). This approximation 
results in diffusive evolution of intermediate waves (e.g., slow magnetosonic, Alfven, and entropy waves) and is the main source 
of dissipation in our simulations. 

The discretization of the computational domain results in cubic computational cells with sides Air = Ay = Az = Al = L/N. 
We adopt standard finite volume index notation to associate variables with the cells in the Cartesian grid: the coordinates of the 
geometric center of a cell with index triplet (i, j, fc) are denoted (Xi, yj , Zk). Finite volume variables centered on the geometric 
center are also assigned the index triplet. Superscripts n and n + 1 denote time states, and the time step At increments time 
from t n to t n+1 . For example, the volume-averaged (in angle brackets) x-component of the velocity in the cell at time t n is 
denoted (u x )™- k . Magnetic field components are centered on the faces of computational cells in the CT method. For example, 
the a;-component of the area-averaged magnetic field, centered on the outer face of cell {i,j, fc) with coordinates (x i+ i , yj , Zk) at 
time t n , is denoted (B x ) n k . , . (For uniform grid spacing we have x i+ i = Xi + Ax/2.) 

An integration of the magnetic induction equation over the cell face with normal parallel to the cc-coordinate direction and 
time interval At results in (after application of Stoke's theorem and replacing time-integrals of electric field components with the 
rectangle rule) the time-explicit finite volume update formula for the area-averaged x-component of the magnetic field 
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and the line-averaged z-component of the electric field (centered on the cell-edge) is 

(E z )? + i tj+hk = ^ / 2 E z (x i+ r,y j+h z,t n )dz. 
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Update formulae for the other magnetic field components are obtained in an analogous manner. The update given by Eq. ( |A1[ | is 
exactly the forward Euler method and results in first-order temporal accuracy. We obtain second-order temporal accuracy with a 
Runge-Kutta method (e.g., |Shu |1997[ >. 

The key to stable and accurate magnetic field evolution with the CT method is the specification of the edge-centered electric 
field components. The z-component of the edge centered electric field with spatial coordinates (x i+ i , y ,- + i , Zk) is computed 
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with an HLL-type formula (cf. Londrillo & Del Zanna 12004} 
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where a± = max(0, ±A±'™ ±A±' NW ,±A±< SE ,±A±< NE ), and A 



± _ 



± c{. are the characteristic wave speeds associated with 



the fast magnetosonic modes. The maximum is taken over wave speeds computed in the four cells sharing the edge indexed 
(i + |,j + |, k), which are denoted with superscripts SW, cell (i,j, k); NW, cell (i,j + 1, k); SE, cell (i + 1, j, k); and NE, cell 
(i + l,j + l,k), respectively. (We have simplified the notation in Eq. ( A4 1 by replacing the subscript indices i+ | , j + I , fc with 
p.) For example, for first-order spatial accuracy we simply set [E% 



*+5,i+,,fe 



= B, 



i+±.j,k\ a yn,j,k 



B 



id,k\- LJ y/i,j+^,k' 



Magnetic field components are centered on cell faces, and cells indexed (i + 1, j, k) and (i + l,j + 1, k) share [By] 



(B„ 



which is assigned superscript E. Similarly, cells (i, j, k) and (i, j + 1, k) share [By] 



B„ 



which is assigned superscript W. To simplify the presentation, we only briefly discuss the first-order scheme. However, all our 
calculations are done with a second-order scheme. For second-order spatial accuracy we use monotonic linear in terpolation (via 
the multivariable minmod limiter) to reconstruct va riab les to the appropriate edges (e.g., Kurganov et al. 2001, and references 
therein), and evaluate the electric field thro ugh Eq. ( A4 1, 

The HLL electric field given by Eq. ( A4 1 contains explicit dissipation due to the second and third terms on the right-hand-side. 
(We refer to the first term on the right-hand side as the ideal part of the electric field.) The dissipation terms mimic the effect of 
physical resistivity due to a non-ideal electric field — u x B + 77J. This becomes evident by considering a subsonic flow with a 
weak magnetic field (appropriate for the turbulent post-shock flows in our simulations). Then, A J « ±c[ « ± cs an d af ~ cs, 



where eg is the sound speed. (Similarly we have a± « cs-) With these approximations the electric field in Eq. ( A4 1 becomes 
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where the "numerical resistivity" 77 mim = | cs Al /i Q decreases linearly with Al 

the parenthesis of the second term on the right-hand side of Eq. ( A5 1 is a numerical approximation for ( V 



In particular, we note that the terms inside 
x B) z = fj, J z at 



(x i+ i , Vj+i , Zk)- (Thus, the analogy to non-ideal MHD with scalar resistivity used in the discussion in Section 



3.7 



is appropriate.) 
However, this 



These terms prevent growth of grid scale oscillations and act to stabilize the evolution of the magnetic field" 

numerical dissipation becomes non-negligible when the magnetic field develops a flux rope structure with the flux rope thickness 

comparable to a few grid cells. 

To investigate the amount of magnetic energy dissipation in our simulations we have initiated a simulation with data from one 
of our models (model B1 2L0 .0E00 at t — 900 ms) and evolved the magnetic field with the induction equation with only the 
dissipative terms in Eq. (A4i (the second and third term on the right-hand-side) retained. The fluid variables (p, u, e in „ etc.) 
are held fixed to their initial values (t = 900 ms). The results are plotted in Figur e [TTl where we plot the magnetic energy E m . i% 
(black lines) and the characteristic spatial scale of the magnetic field A mi ,„ (cf. Eq. (fTTT and proceeding text; gray curves) versus 
time. We see that the magnetic energy decays rapidly with time in the absence of theldeal part of the electric field (black dashed 
curve). Initially, it decays on a millisecond time scale (the decay rate is about 380 s _1 at t = 900 ms). The characteristic spatial 
scale of the magnetic field increases as a resuh of diffusion (gray dashed curve), which results in a decrease in the decay rate. 
The magnetic energy increases with time, and A mag stays roughly constant in the full run (solid lines). These results confirm our 



claim in Section |3T6| that numerical diffusion plays an important role in our simulations, and support our claims in Section 3.7 



that our simulations grossly underestimate the magnetic energy growth rate due to turbulence-driven magnetic field amplification 
in CCSNe as a result of the SASI. Thus, we expect the magnetic energy to grow on millisecond time scales, if dissipative effects 
can be ignored in the supernova environment. 

The diffusive evolution of intermediate waves by the HLL Riemann solver motivates us to improve our MHD scheme by 
including more w ave families in the Riemann solver in future applications. In particular, the HLLD Riemann solver (Miyoshi 
& Kusano~||2005 ) includes Al fven and entropy mod es and is an attractive option (perhaps in combination with an improved 
interpolation scheme). Indeed, Mignon e et al. | (2009) have developed an HLLD-type scheme for relativistic MHD, which shows 
significantly improved resolution on small scales when compared with the corresponding HLL-type scheme. 
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Figure 17. Evolution of the total magnetic energy (black lines) and the characteristic spatial scale of the magnetic field A mag (gray lines; cf. Eq. jllf and 
proceeding text). Solid lines represent data from model Bf 2L0.0E00, while dashed lines are from a simulation initiated from model B12L0.0E00 at t = yuO ms 
and evolved with the i nduc tion equation for ten milliseconds, with the fluid held fixed to the state at t = 900 ms, and only the second and third terms on the 
right-hand side of Eq. (A4l included. 



